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ABSTRACT 


Large Reynolds number approximations are examined 
for certain internal separated flows by considering -various 
limits. It is argued that the equation governing the flows 
except possibly for some subregions, has the form of the 
boundary-layer equation, as the inertial and viscous forces 
are of the same order. A method of calculation is then 
developed which uses the eigenfunctions of the linear 
problem governing the far-downstream flow. It leads to a 
set of coupled, first order, ordinary differential equations 
It is first tested for two cases for which results based on 
the Navier-Stokes equations are available. Details of the 
flow including streamlines in the recirculating regions are 
well predicted. It is concluded that the equation represent 
the large Reynolds number behaviour of the flows considered. 
Flow problems with three types of wall geometry are then 
solved. They are a channel having symmetric sudden expansio 
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a channel having asymmetric sudden expansion and a "base 
placed symmetrically in a channel. Calculations are 
performed for various values of the ratio of the upstream 
channel width to that downstream. Also two types of entry 
conditions representing a long and a short inlet are consi- 
dered. Results which are presented include centreline 
velocity, velocity profiles, development length, streamlines, 
skin friction, pressure distribution and pressure recovery 


coefficient . 



CHAPTER i 


INTRODUCTION 

Separated flows constitute an important tut difficult 
area of fluid mechanics. Bound ary- layer approximations 
provide substantial simplifications in attached flows at 
large Reynolds number. In particular, the outer flow can be 
obtained to the lowest order, by an inviscid calculation in 
which the condition of zero normal velocity at the wall is 
satisfied. This is because the bound ary- layer and the wake 
where the viscous forces and the inertial forces are of the 
same order, are' rather thin and their effect on the outer 
flow is of a higher order. Separated flows, on the other 
hand , have quite thick rotational regions and the outer flow 
is affected substantially. Consequently, the outer flow 
calculation presents certain difficulties which have yet to 
be fully resolved. 

The present thesis deals with a simpler clans of 
separated flows which has one ingredient in common with 
external separated flows, namely that they have thick 
rotational regions. We consider the flows confined between 
two walls -(figure 1), which belong to the category of 
internal separated flows. One such flow is the flow in a 
channel having a sudden expansion in which recirculating 
regions are formed downstream of the expansion. As we shall 
see in Chapter 2, there are two types of internal separated 
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flows which, differ in their streams? ise length scales. 


limiting equations for various longitudinal and 
transverse length scales are examined as B-^no (Cnapter 2). 

It turns out that for the type of internal separated flow 
problem referred to earlier, the streams? ise length scale is 
of order E. As the inertial forces are of the order of the 
viscous forces in the limit considered, the equation has the 
form of the -hound ary-layer equation. However, the longi- 
tudinal and transverse length scales of the present problem 
are different from the classical boundary-layer scales. 

Also, the pressure gradient is an additional unknown which 
is determined by using four boundary conditions arising from 
the walls unlike the boundary-layer problem. 

A calculation procedure of the equation, which is 
developed in Chapter 3, exploits the ultimate development of 
the flow into Poiseuille flow by taking the solution in 

the form 


vU 


(x ,y;R) = v U(y) + > a. (x/E) (y) 

■■■— 1 I j_ 

1 


where x and y are the streamwise and transverse coordinates 
and is the stream function. Usually in an integral 
relations method , the expansion set and the weighting 
functions are chosen for their analytical properties and 
computational simplicity. “e however, adopt .a. different 
approach here. The flow at large downstream distances is 
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a small perturbation of the Poiseuille flow which leads to 
an eigenvalue problem. Its eigenfunctions form an efficient 
expansion set as the resulting set of coupled ordinary 
differential e! ua tions tend to get decoupled with increasing 
x. The -weighting functions are the natural ones arising 
from the eigenvalue problem. The given flow problem then 
reduces to integrating these differential equations which 
are coupled quasi-linear and of first order. The initial 
condition has to be obtained by matching the far -flow field 
with the ne ar—f field. In the absence of any rigorous 
description of the near-flow field, a hypothesis is made to 
give the initial condition. 

The sche 1116 ° f calculation is first applied to two test 
cases for which numerical results based on the complete 
lavier-Stoke s equations are available, and can be compared 
with the present calculation. The first problem to be taker 
up is the now classical entry flow problem (Chapter 4-). Its 
simple geometry forms a suitable testing ground for theories 
and methods and hence Is considered here, although no 
recirculating regions are present. The second test case is 
the flow in a. channel having a symmetric sudden expansion 
with an expansion ratio of two. Comparison of the present 
calculation wi'th that based on the Navier-Stoke s equations 
for large R (Hnng and Macagno , 1966) is found to be very ' 

good . Not onXy the gross quantities of interest such as thf 
length of the recirculating region and the total recirculst; 
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are found to agree well, "but also the streamline pattern is 
•well predicted. These results are discussed in Chapter 5. 

The above comparison shows that the present method of 
expansion into eigenfunctions of the far downstream flow 
provides a suitable calculation procedure , and it also 
demonstrates the applicability of the limiting equation. 

We thus conclude that the streamwise length scale for such 
flows is of order R, the viscous term and the inertia term 
are of the same order, and the governing equation is 
parabolic despite the presence of recirculating regions. 

Flow problems with three different types of wall 
geometry are solved by using the present method (see 
Figure 1). One type illustrates the applicability of the 
method to asymmetric flow problems. Another type models a 
base between two walls. Calculations are made for various 
values of the geometrical parameter h, and are included in 
Chapter 5. The results of these calculations are believed 
to be new. 

As the method uses an eigenfunction expansion, the 
calculation, in principle, should be continued with increas- 
ing number of terms until a satisfactory convergent trend is 
obtained. However, we have used a rather small number of 
terms in the calculations. For example, the favourable 
comparison in the ca.se of the symmetric sudden 



expansion mentioned earlier is obtained with only three 
terms. The maximum number of terms used in the present 
calculations is therefore limited to five. Obviously the 
choice of the number of terms would depend on the wall 
geometry, the entry profile and the desired accuracy. 

If the number of terms to be used in a calculation is 
fixed, a difficulty may be encountered, in the integration 
of the ordinary differential equations. This difficulty is 
essentially due to the non-linear evolution term and is 
a general feature whenever a moment method or an integral 
relations method is applied to the present type of partial 
differential equation. This issue is discussed in Section 
3.5. Phase plane analysis of the ord inery d ifferent ial 
equation when two terms are used , is given to illustrate 
this feature. 

We have restricted our attention to steady laminar 
two-dimensional flow of a Newtonian incompressible fluid, 
although the basic elements of the analysis would be appli- 
cable to a larger class of problems. We have however chosen 
to apply the method to three types of wall geometry, for a 
range of the parameter h, two type of entry conditions, and 
three choices of the number of terms with a view to obtain 
insight into the performance and the applicability of the 
method . 



CHAPTER 2 


LIMIT EQUATION FOR LARGE REYNOLDS NUMBER 

2.1 GOVERNING EQUATION 

The equation governing the two-dimensional plane 
steady laminar motion of a Newtonian incompressible fluid 
is 

j( y 2 y> ) = r _1 v 2 \7 2 , (2-1) 

where ^ is the stream function. The two cartesian coord i- 

2 

nates are denoted by x ' and y, and \7 and J denote the 
Laplacian and the Jacobian operators respectively. The 
above equation is written in the non-dimensional form, R 
being the Reynolds number. Equation (2.1), known as the 
vorticity transport equation, is obtained from the two 
momentum equations and the equation of continuity. The 
above equation, without any approximation , will be denoted 
by NSE . 

2.2 LARGE REYNOLDS NUMBER LIMITS 

Let the scales of x, y and be denoted by 1 , 1 

x y* 

and / respectively. Then the orders of the various terms 
in equation (2.1) are 

( //b^Ly, ^ 2 /L x L^), ( Y/RL x , (f/RL 2 ! 2 , ^/RLy) . 

The first bracket shows the orders of inertia, terms, while 
the orders of viscous terms are shown in the second bracket. 
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Terms of the lowest order in either group will depend on 

the order relation between the two length scales I and L . 

x j 

Therefore three sets of limits are possible depending on 

whether I v ,'j L r , L c_ I- , or I c: I . Possible limit 
x ^ y x x j 

equations are given below and each limit equation has been 
designated by E, V, B or IT standing for inviscid , viscous, 
boundary-layer or Navier -Stokes types of limits. 


E 


IT 


Y 


E 

y 


B 

y 


J(v 2 'f , y ) = 0 

J(\7 2 'V , S ) = R _1 Y? 2 \7 2 'V 

} > 

p p 

x7 ^7^=0 


1 C. 1 

y x 


j( S' , o ) = o 

■ yy 


j( vjp S ) 

< yy 1 


P 1 


’ yyyy 



KjJ 

1 yyyy 


o 


The order symbols used here have the following meaning. 

1 n^_ i 0 <Ilim l /jj oQ 

x * e-co x y 

W ^ K * lim 1/1 = 0 

3' a-^-c,- x I 


x 


y 
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I C 1 
X "V 


j( u> , t> ) = o 

V i XX i 

J( V , ) = R _1 ib ^ 

v t xx 5 J ! xxxx 


xxxx 


( 2 . 2 ) 


Suffixes x and y following ^ denote respective partial 
derivat ives . 


limit equations E end Y can "be obtained from B , 

7 1 7 

and E and Y from B . Limit equations E and E are 
x x x y x 

contained in E and V contains Y^ and Y . Distinguished 
limits are those corresponding to E, Y, B and B and E is 


the complete equation of motion. 


E is the Euler equation governing inviscid fluid flow. 
Stokes equation in which only viscous terms are retained is 
designated by Y. Interestingly, Stokes equation is obtained 
for large Reynolds number E. However, the appropriate 
Reynolds number i(R goes to zero. B„ t is the familiar 
bound ary- layer equation differentiated once with respect to 
7 to eliminate the pressure term. E is sometimes referred 
to as the inviscid bound ary- layer equation. B and E_ are 

X X 

x counterparts of B and E respectiveljr. 


The domains of validity of the various limit equations 
can be readily visualised with the help of a three-dimen- 
sional plot with axes 1, , 1 and y (figure 2). This plot 

x y 
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is of a special nature as the axes are order axes. 'That is, 

a point, say 1), denotes L rj R, I ■ - VE and V -v» 1. 

x y 

The choice of a log scale for the axes makes the represen- 
tation simpler. A suitable unit for the plot is R. In 
.Figure 2, OA = OB = OP = R. 


1 L represents the plane IGHJ in Figure 2. 1 l/R 
x ' y 

on the line EF in this plane, and the limit equation IT is 
■valid on this line. The domains of validity of the limit 
equations E and V are the portion above and below the line 
EF of the plane L„ fu that is EGHF and EEFJ respectively, 
excluding the line EF in both cases. 


As equation (2.1) is symmetric in x and y, the order 

space plot is symmetric about the plane 1 ru I • Consider 

x y 

the case when 1 0 L . The order of the inertia term is 
y x 

2 3 i 

1 /Il_ 1 and that of the viscous term is 1 /HI. When 

X J y 

/ — 1 

the two terms are of the same order , y 1 /L ru B. which 

represents the plane EKLF due to the choice of log scales 

for the axes. This plane, to be called B^_, meets the plane 

J 

1 nj 1 along the line EF. Since 1 c_ L , B is restricted 
x y y x y 

to only one side of the plane I > o L . limit eauation B 

x y ^ y 

is valid on this plane except for the line EF. The domain 
of validity of the limit equation E and Y. are- regions 


above and below the plane B. such that L 


that is. 


regions between the planes BGHF and EK1E and between the 
planes EIJE and EKIF respectively. A similar construction 
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is carried out when 1 C, 1 . Each plane in the figure is 

x y 

infinite or semi-infinite. The construction of the figure 
is restricted tolSbl for clarity. 

The line joining 0 and Q in Figure 2 is in the plane 
I ^ 1 and such that u = w ru "tf/L 1 . OQ meets the line 

x i y y 

of intersection FC of the plane L.m 1 and the plane 

IK. 

B r (EKLF) , at Q. Therefore the point Q represents the 
classical hound ary- layer limit. Order coordinates of the 

i i 

point Q are (1, R 2 , R 2 ). 

It is clearly seen from Figure 2 which shows the 

possible limits of HSE, that although there are infinitely 

many limit processes (each point of the figure corresponds 

to a limit process) the possible outcomes are finite in 

number. The domain of validity of the limit equations E, Y, 

B and B are planes whereas those for E , E , V and Y are 

three-dimensional regions. Each of these regions is bounded 

by two of the four planes BGHF, EKIF, EIJF, and EMUF. For 

example, the domain of validity of V is the region between 

the planes EMF and EIJF. Therefore limit equations E^, E , 

V and Y are each contained in the limit equations corres- 
x y 

ponding to the planes which bound their respective domains 

of validity. That is Y is contained in Y and B and so on. 

x x 

The four planes mentioned above meet along the line EF 
which corresponds to USE. The limit equations corresponding 
to these planes are therefore principal limit equations of 
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NSE. B and B. differ only in the definition of the 

^ y 

coordinate axes. Thus we have only three principal limit - 
equations of USE. They are the Euler equation, the Stokes 
equation and a third equation with which we associate the 
classical hound ary- laye r . However, this limit equation is 
more general. 

The above construction of the order space helps in 
visualising various limits of USE. A somewhat different 
construction of the order space can he one where we choose 
the order axes to be 1 , L and instead of V the order of 

x y 

the streamwise velocity. This, however, represents a 
linear transformation in the order space. 

2.3 PHYSICAL PROBLEM 

The flow is considered to be confined between two 
walls. A simple example of an internal separated flow is 
the flow in a channel having a symmetric sudden expansion 
(Figure 1). Recirculating regions are formed downstream of 
the expansion, and flow reattaches some distance downstream 
and finally evolves into a. Poiseuille flow. An essential 
feature of this problem is that the flow is bounded in the 
transverse direction and, therefore, this length scale can 
be taken to be of order one. Two types of internal 
separated flows are identified later in Section 2.6. 

In external separated flow, the interplay of the 
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recirculating flow and the outer irrotational flow requires 
that the shape of the "bubble or the "boundary of the 
irrotational flow be determined as a part of the problem. 

As a result the transverse length scale of the separated 
wake is difficult to ascertain. However, the boundary 
conditions on the two walls are sufficient to formulate the 
present pr o ble m . 

2.4 LIMIT EQUATIONS WHEN L C O 1 

3 ^ 

When L 1, (i.e. the plane normal to L axis in 

y y 

Figure 2) possible limit equations with their domain of 

validity are shown in Figure 3. The limit equation N is 

valid at P. E and V are valid on the axis above and below 

P respectively. B is valid on PA whereas B on PD. Domains 
J y x 

of validity of E E . V and V are regions OPA, OPD, APT 

j x y x 

and DPT respect ively. line OG in Figure 3 is such that 
v = - rv- 7// L x 1. Also, along the 1 axis u .•■*/ 1. There- 
fore, if the flow velocities are bounded for large R, the 
region of interest is that below AOG in Figure 3. Point G 
corresponds to the usual boundary-layer when the v axis is 
along the surface . 

For the physical problems considered the total flow 
across any section remains constant and in its non-dimen- 
sional form is independent of R. Therefore , J ,v 1. From 
Figure 3 we see that there are then two principal limits. 
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They correspond to point 0 and point A. 

The limit equation corresponding to point 0(1, 1,1) 
is the Euler equation 

J ( V > V ) = 0 (2.3) 

The flow in this limit is inviscid . The above equation is 
elliptic in nature. Consequently, if the solution is 
disturbed at any point, its effect is felt at all other 
points of the region. 


The limit equation correspond ing to point A (R ,1 ,1) 
in Eigure 3 is 


\U 

>y 


r yyx 


O' U.-' 

'X iyyy 



■ yyyy 


(2.4) 


The length scale in the x -direction is of order R. The 
above equation is parabolic so that there is a direction at 
every point, called the characteristic direction. The 
effect of a perturbation introduced at any point in the 
region is sensed along this loca.1 direction. 

2.5 DEVEIO MEN T 0E POISBUILIE EOT 


We first consider the development ox Eoiseuille flow 
as a prelude to the internal separated flow problems, as 
the flow between the two parallel walls develops ultimately 
into Eoiseuille flow. Wilson (1969) has considered the 
development of Eoiseuille flow using NSE. A small 
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perturbation analysis is carried out with the perturbation 
ass ume d in the form <i*(y) exp(-Ax), where x and y are the 
streamwise and transverse coordinates. An eigenvalue 
■problem is obtained for -£> , X being the eigenvalue. For 
large E, two sequences of eigenvalues are obtained. The 
members of the first sequence arc of order one, while those 
of the second are of order l/E, An eigenvalue of order one 
indicates that the perturbation is effective over a distance 
of order one. The eigenvalues of this class are positive 
as well as negative. The perturbations are, therefore, 
sensed in both the upstream and the downstream directions 
over a distance of order one. The streamwise extent over 
which the perturbations corresponding to the second sequence 
are felt is of order E. These eigenvalues are such that 
the perturbations decay in only the downstream direction. 

It then follows that if a small disturbance is 
introduced on the Poiseuille flow, it propogates both 
upstream as well as downstream. However , the extent of 
upstream propagation is of order one and that downstream is 
of order E. Also, it is the second mechanism which is 
dominant in the downstream propagation of any disturbance as 
its decay rate is considerably slower. 

The eigenvalue problem for the second type of eigen- 
function can be derived from the limit equation (2.4). 
Qualitatively the mechanism -of propagation of the 
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disturbance should therefore "be similar to that of equation 
(2.4). Since the equation trs a, characteristic direction, 
the disturbances are expected to propagate in only one 
d irection. 

It can then be argued on the basis of the nature and 
length scales of limit equations, and the analysis of Wilson 
that there are two modes by which disturbances are propa- 
gated through the flow. If there is an upstream influence 
then this flow is expected to be governed by the limit 
equation (2.3). The second mechanism is dominant for down- 
stream flow, and is described bj* the limit equation (2.4). 

2.6 TWO TYPES OF EITEE1TA1 SIIilUTZD FIOWS 

Two simple geometrical arrangements causing separation 
in internal flows are sudden contraction and sudden expan- 
sion (see Figure 4) . The separated flows in these configu- 
rations are qualitatively d if ferent . In the case of sudden 
contraction, the flow separates upstream of the geometrical 
configuration responsible for it. The streamwise length 
scale is, therefore, of order one, and the separation bubble 
appearing ahead of the contraction will be of order one. 

The second type of separated flow is that downstream 
of a sudden expansion. The dominant mode of propagation of 
the effect of the sudden expansion in the downstream 
direction is that corresponding to equation (2.4). The 
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streamwise length scale is of order E, "that is a typical 
length; say, the length of the separated bubble, behaves as 
r. Viscous diffusion over a distance R is significant . 
Therefore, viscous forces are not negligible in the dynamics 
of the flow in the bubble. The appropriate limit equation 
is (2.4). Separated flows of this class are dealt, with in 
t he or t.- se nt w 0 r k . 


The above rather intuitive conclusion can also be 
reached by the following argument. When 1 ru 1 end L '"- y 1, 

x y 

viscous terms are of higher order. In the 'near 1 limit, the 
flow is governed by the Euler equation. Discontinuous 
surfaces ere allowed. One such surface is along the sepa- 
rating streamline. A vorticity layer is required for the 

vorticity to diffuse in the bubble. The thickness of this 

l/3 

layer is of order (x/R) and therefore for any fixed x, 
it shrinks to zero in the limit. When x R, the vorticity 
layer thickness is of the order of the transverse length 
scale of the channel itself. The ’near’ limit is, therefore, 
inapplicable when x ^ R. A ’far 1 limit given by ^ R is 
required which leads to equation (2.4). 

2.7 EARLIER RESULTS OR STREAMWISE LENGTH SCALE 

We examine the behaviour of a typical streamwise 
distance for large Reynolds number in the earlier numerical 
calculations with a view to test the above conclusion 
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regarding the nature of separated flow. 


Hung end Macagno (1966) have calculated internal 
separated flow in a. channel having a symmetric sudden ex- 
pansion. Their calculation is "based on i'TSB and uses a 
finite -difference scheme. Two methods have been used to 
study the problem. In the first, steads 7 state equations 
are solved, while the se cor’d approach is based on unsteady 
equations and the calculation is continued till a steady 
state is reached. The second approach is considered to have 
advantages ever the first for large R. The ratio of expan- 
sion (downstream to upstream width) is two and the calcula- 
tion is performed for E (based on the average velocity and 
the half-width of the channel downstream of the expansion) 
upto about 160. Flow slightly upstream of the expansion is 
taken to be fully developed i.e. a Poiseuille flow. Two 


typical streamwise lengths are shown as functions of R in 
Figure 5. Thej r are, the distance of the point of re attach- 
ment and the distance of the centre of the ’ eddy’ measured 
from the expansion. They clearly show a linear trend for 
sufficiently large R. This problem is also considered by 
Morihara (1972) . Calculations of RSE are performed upto 
R = 20, using a different numerical scheme. The length of 
the 'eddy' obtained by Morihara is also shown in the figure . 
Kis results also show a linear trend. Flow in a. gradually 
widening channel has been calculated by Dorodnitsyn and 
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Me Her (1970). 'Their Figure 6 
reversed flow region which increases 
linearly with R. 


approx imately 


We thus conclude that the earlier calculations based 
on iTSE support the view that the stre atnwise length scale is 
of order R and it is then expected that equation (2.4-) 
should govern the large Reynolds number flow except possibly 
for some small sub-regions. 


2.8 FORMULATION OF THE 1R0BLEM 


We shall consider the flow in a sudden expansion as a 
prototype lor our discussion, although the conclusions are 
expected to be more general and applicable to other similar 
problems such as shown in Figure 1. 

The flow downstream of the expansion is divided into 
two sub-regions corresponding to 1., •• v* 1 and L oj R. The 
limit equation for the ’far’ region (L„ 'v/ R) is (2.4). As 

•A. 

a result of the limiting process, three derivatives in x are 
lost. Equation (2.4) can therefore satisfy only one condi- 
tion in x which is obtained on matching with the solution for 
the ’near' region (l^n. 1). The ’near’ region solution is 

u\. 

somewhat difficult to construct. As the limit equation for 
this region is elliptic in nature , it would further require 
matching with the far upstream flow. To obtain the initial 
condition for equation (2.4) in the absence of matching , we 
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make a working hypothesis which is to be tested subsequently, 
that the initial velocity profile is the same as the far 
upstream flow. 


We now eliminate R in equation (2.4) by scaling the 
x coordinate b j R. It then becomes 




■ y 


yyx 


■4/ Ur 

■x iyyy 


JJJJ 


(2.5) 


Unless otherwise stated , x would henceforth imply the 
contracted streamwise coordinate. In this work R is taken 
to be based on the average velocity and the half -width of 
the channel downstream of the expansion. 


The initial condition required for the integration of 
equation (2.5) is given as 

tv-(x = 0, y) = ( y ) . (2.6) 

This condition is equivalent to prescribing the streamwise 
velocity at the initial section. The no-slip condition at 
the walls requires 

r'= ± 1 > U; r = 0 at y = ± 1 . (2.7) 

In the coordinate system x (contracted) end y, the 
'near 1 region shrinks to a plane, x = 0. Equation (2.5) ear- 
then be regarded as the field equation valid for x> 0. ft 
x = 0 an initial condition (2.6) is prescribed. Equation 
(2.5) along with conditions (2.6) and ( 2 . 7 ) constitute 
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the mathematical problem and subsequently we test if they 
form a. sufficient basis to calculate the flow. 

Equation (2.5) can be integrated once with respect to 
y. The constant of integration is the pressure gradient. 
However, as we will see, it is preferable to work with the 
higher order equation (2.5) than with its integrated version. 

2.9 FURTHER COMMENTS 


The equation obtained for the calculation of flow 
development is similar to Prandtl’s bound ary- layer equation 
BEE. However, it should be noted that the viscous region 
has a thickness of order one whereas the thickness of the 

i 

classical bound ary- layer is of order R~ 2 . Also, instead of 
pressure being impressed by the external flow, it has to be 
obtained as a part of the calculation here. However, there 


are four boundary conditions in y to be satisfied. 


.tie 


confinement of the flow can thus be said to generate its own 
pressure field. Strong interaction of the separated region 
with the ma inflow which is a feature of the external 
separated flow is absent in the present internal separated 
flow . 


It may appear surprising that the equation which we 
intend to use for flows having recirculating regions is of 
the boundary-layer type. BEE are known to possess a singu- 
larity at the point of sere skin friction (Gold ste in , •' 1 948) , 
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and it is "believed that the continuation of the calculation 
beyond this point is generally not possible (Brown & 

Stew art son, 1969). Flow reversal occurs beyond the sepa- 
ration point. Since the argument used in deriving BEE 
does not require the absence of reversed flow, BEE permit 
reversed flow as long as the bound ary- layer thickness is 
of order R - ® . Briley (1971) obtained the NS solution for 
the flow in a separation bubble by means of a finite- 
difference scheme, which with increasing R appears to 
approach the corresponding solution of BEE. 

Consider Figure 6 (taken from Brown &. Stewart son, 
1969). The conditions on an initial section 13, the wall 
k F and the edge of the boundary-layer BE are given. Since 
the flow in C-EF is reversed, the convective propagation of 
a disturbance in the sense of decreasing x is possible. 

In such circumstances, the parabolic equation is said to 
be of the mixed type . There seems to be an intimate 
connection between the theory of parabolic equations of the 
mixed type and separated flow. Apparent ly , the requirements 
on the boundary condition are severe for such equations. 

The theory of such equations has not yet been sufficiently 
developed to give conditions which the given data must 
satisfy in order that the problem is well posed and a 
solution exists. In the absence of any formal theory, 
qualitative arguments are advanced for the appropriate 
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■boundary conditions. It may appear tnat in addition bo the 
already specified conditions on AB, BD and AE, a condition 
on EE (Bigure 6) is required in order to determine the 
solution over ABDB. However , if the boundary conditions 
were to be completely arbitrary, flow upto point G will be 
affected by changing the condition over GB and EE, and 
possibly the separation point G itself will shift" but the 
flow in ABCG uniquely determines the point G. Perhaps some 
condition on EB is needed in general to determine the eigen- 
solutions if any. ‘The above argument also suggests that the 
boundary condition to be prescribed should be in a certain- 
class for the solution to be regular. The regularity of the 
solution imposes certain restrictions on the boundary condi- 
tion, even when the region of interest is ABCG. when the 
given condition is otherwise, the flow manifests a singu- 
lar it y . 

Based on the above discussion, the appearance of a 
singularity can be viewed from two angles. If we are 
looking for 'thin' separation (the contention here is that 
it is regular) , the singularity represents the incompatible 
nature of the supplied data. Cat her all and Mangier (1966) 
designed a thoughtful scheme in which the supplied data had 
partial freedom to adjust in order that the solution be 
continued in a regular waj at each step. On the other hand , 
the singularity would mean the failure of the boundary-layer 
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limit as an approximation to tjpe actual physical phenomenon 
('thick 1 separation). 

There are two distinct issues in the calculation of 
separated flow, namely, the equation and the boundary 
conditions required for the formulation of the mathematical 
problem, and a method of calculation for the formulated 
mathematical problem. Equation (2.5) is considered to be 
appropriate for the present class of problems. The 
boundary conditions prescribed on the walls (2.7) seem to 
be of the simplest type. If any difficulty arises, it can 
be attributed to the initial condition (2.6). ¥e shall 
regard the problem as well posed and assume that a solution 
exists. 


The usual methods of integration of parabolic equa- 
tion are not expected to be applicable in the presence of 
reversed flow. An iterative technique is used by Klemp and 
Acrivos (1972) in which the shape of the reversed flow 
region is assumed , and then the solution in the forward 
direction and backward direction are obtained separately by 
a marching technique suitable for parabolic equations. The 
shape of the reversed flow region is iterated to match the 
flow quantities across it. In solving a parabolic equation 
of mixed type, Ban and Euerti (1969) have used a numerical 
scheme which is normally applied to elliptic equations. 
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Integral relations methods or methods of weighted 
residuals are well known for solving partial differential 
equations. Their convergence can he shown under certain 
circumstances. The method of calculation used in the 
present work to solve equations (2.5 - 7) belongs to this 
class and is developed in the next chapter. However, the 
approximating and weighting functions are eigenfunct ions 
of a related problem, unlike the usual functions which are 
often selected for their analytical properties and 
computational simplicity. 



CHAPTER 3 

METHOD OP CALCULATION 

3.1 INTRODUCTION 

This chapter deals 'with the formulation of a method 
to calculate the flow development "between two parallel walls, 
which can "be applied to internal separated flows. The 
governing equation and the required "boundary conditions are 
given in Section 2.8. 

Since the flow evolves ultimately into Poiseuille 
flow, equation (2.5) can "be linearised about this flow for 
sufficiently large downstream distances. It then leads to 
an eigenvalue problem and the perturbations decay exponen- 
tially with increasing streamwise distance. In the region 
where the flow differs markedly from Poiseuille flow, the 
linear approximation does not hold , and nonlinear terms must 
be retained. However, if the nonlinear interaction terms 
between the modes decay rapidly, an expansion in the eigen- 
functions could lead to an efficient computing procedure. 

3.2 LINEAR PROBLEM 

As x-^co , the channel flow develops into Poiseuille 


flow given by 
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Let the flow he expressed as a sum of Poiseuille 
flow and a perturbation , given by 


Vu + 

>p 1 1 


(3.2) 


Substituting (3.2) in equation (2.5), we get 

kb 

Tyyyy >p nyyx 


uv vi/ + V w k H'ly l Hyyx ~ Tlx ''H yyv * ^ ^ 
i p slyyx 1 p 'lx J J J - 


Lash indicates differentiation with respect to the argument. 
Neglecting the quadratic terms in Vi/ and substituting for 
ypp from (3.1), we get the linearised equation 

Tiyyyy ” 5//2 ^ “ y ) 'i"iyyx 3 Hix = 0 * (5*4) 

The form of the above equation suggests a solution of the 
form 

= <P( y) e~ Xx , (3.5) 

for some X . Substituting (3.5) in (3.4) we get 

4p 1V + X [3/2(1 - y 2 ) 4" + 3=p] = 0 , (3.6a) 

where X is an eigenvalue. The wall conditions require that 

4> = 4=*’ = 0 at y = + 1 . (3.6b) 

Linearization about the Poiseuille flow thus leads to an 
eigenvalue problem EVP, and the perturbations decay expo- 
nentially in the positive x direction, if X is positive. 
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3.3 THE EIGEF7ALUE PROBLEM 

The form of (3.6) suggests that its solutions can he 
separated into odd and even solutions. An eigenfunction is 
called odd or even depending on the antisymmetry or the 
symmetry of the associated streamwise velocity component. 

EVP (3.6) has been studied earlier by Gillis and 
Brandt (1964) and by Wilson (1969). Gillis and Brandt have 
made a direct numerical computation of even eigenfunctions. 
Wilson has obtained equation (3.6) as the equation governing 
the leading term in a large Reynolds number expansion. 

Equation (3.6) can be integrated once to give 

b + X [3/2(1 - y 2 ) O + 3y<t>|= Constant . (3.7) 

The constant of integration in (3.7), when multiplied by 
e~ gives the associated pressure gradient. Evaluating 
the left side of (3.7) at y = 0, the constant of integra- 
tion turns out to be zero for odd eigenfunctions. Therefore 
the pressure field associated with an odd eigenfunction is 
uniform. 

It has been shown by Wilson that the odd eigenvalues 
are real. It has not been possible to prove a similar 
result for even eigenvalues. However all the eigenvalues 
calculated by using numerical methods are real (Wilson, 1969). 
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The real eigenvalues of (3.6) can he shown to he positive 


hy multiplying (3.6) hy <p and 
leads to 

^r- 1 

I dy - 3/2 a j 1(1 - 

U ! --1 


integrating hy parts, which 


2 ,2 
y ) dy 


1 

' , 2 . 
j ay 

-1 


= 0. (3.8) 


Now consider 

1 ,2 1 

X = 1(1 - y 2 ) 4 s ' dy/ U> 2 dy 

A ' ^ 

Stationary value of ;/ leads to Legendre's equation. The 
minimum value of ^ , therefore, is two. Hence the term in 
the bracket of (3.8) is positive. Therefore X is positive. 


Wilson has reported that machine calculations have 
shown no real eigenvalues in the range -200 < X < 0. 
Although it has not been definitively established that all 
the eigenvalues are real, it is conjectured that they are 
real and therefore positive. 

Appendix 1 contains an outline of a scheme for 
calculation of eigenvalues and eigenfunctions . Fourth order 
Runge-Kutta-Gills method has been used to integrate the 
differential equation (3.6). Eigenfunctions of higher order 
oscillate more rapidly and consequently have larger deriva- 
tives, which require a suitable method of integration. 
Estimates of A for higher modes have been obtained in 
Section 3.3.1. Table 1 contains the first few eigenvalues. 
They are numbered in increasing order and are alternatively 
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odd and even. The lowest eigenvalue is odd in nature, and 
hence the velocity perturbation which persists the longest 
is anti symmetrical. 


m 

1 

2 

3 

4 

5 


m 

14.4535 

18.8146 

48.8857 

57.5463 

104.572 


TABLE 1 EIGENVALUES 0? EQUATION (3.6) 


3.3.1 ESTIMATES EOR LARGE Ox 

Analysis of equation (3.6) for large \ is somewhat 

similar to that of the Orr-Sommerfeld equation. There are 

four independent solutions, two of which correspond to 

length scale one, and are denoted by n and b . The 

'A ^ "B 

differential equation governing cp and ci> , in the 

■ l A "B 

limit A — > cc , is 

(1 - y 2 ) <£>*' + 2cm = 0 . (3.9) 

One of the solutions, , is a regular function. However 

singular at the wall (log singularity) and therefore 
an appropriate ’wall solution’ is required. 
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The other two solutions are rapidly varying . These 
solutions are obtained using the WKBJ method . he denote 
these solutions by and . They are obtained by 

writing 


^'n d = exp ^.' gdy ^ * 

We assume a series solution of g in powers of a 2 


( 3 . 10 ) 


g 


_1 



+ g. 


+ 


(3.11) 


The series solution so obtained is 

— 

g = + i >? /3?2 (1-y 2 ) 2 + 5/2 y(1-y 2 ) + ... , (3.12) 

which is valid only if 


U /3 (i - y 2 ) I ?> 1 


(3.13) 


Substituting (3.12) in (3.10)? we get 


2 \ ^ ^ ^ Cos \ ± f n. /<~> \ \ \ t II A ,-2 \ , (?• ^ f -z. A A \ 


d-y‘ -) Sin |K3/2>0 2 (y /(1-y ) + sin y) ! .(3.14) 


Condition (3.13) is violated near the walls and therefore 
the above solutions are invalid there. The proper stretch- 
ing for the solutions near the wall can be obtained by 
considering the limit of (3.14) when y approaches the wall, 
let 


s = 1 + y 


(3.15) 


Then as s -*• 0, the limit of (3.14) is 
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4^— ( 2 s )" 5/4 ^[-|(3/2/\)^ (-^T/2 + 4/2/3. s 3/2 )j . (3.16) 

The proper stretching required near the wall, therefore, is 
-1 /3 

X . This can also he obtained by considering equation 
(3.6) or (3.13). The 'near wall' solutions are obtained in 
terms of Airy's function and matched with the WKBJ solutions, 
which gives 

\mn 2 . (3.17) 

Detailed calculation leads to 

Odd ^ X n ^32/3 (n + 1/12) 2 , (3.18a) 

Even s > n ~ 32/3 (n + 1/4) 2 . (3.18b) 

It should be noted that the leading term of eigen- 
values are real and positive, and the eigenvalues are 
countably infinite. Dor a given n, the odd eigenvalue is 
smaller than the even. If the eigenvalues are arranged in 
the increasing order, they are alternately odd and even. 


3.3.2 ORTHOGONALITY RELATIONS 


The present eigenvalue problem is of the form 


ci>+ A jV\ db = 0 , 


,J ~ -fn _ n 


(3.19) 

( 3 . 20 ) 


. . . 

and <J u 'u are the linear operator 


where 
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X = dV<3y^ , = 3/2(1 - y 2 ) d 2 /dy 2 + 3 , (3.21a) 


and X represents the boundary conditions 


= 0 ; 4‘ = -i 5 = 0 at y = + 1 


( 3.21b) 


The eigenvalue problem encountered here is not of the usual 
type. It can be considered to be of a more general nature 
and reduces to an ordinary eigenvalue problem, whenjVb is 
one. We shall assume that the eigenvalues are real and 
therefore positive and denumerably infinite. Now a scalar 
product is defined by 


1 

(<XG ) = i 

-1 


( 3 . 22 ) 


Corresponding to X and ,JVb , we define the adjoint operators 
by the relation 


8 Z>C> - 4>X*0 = dP^AbW )/dy , 

'H - rp>: *« = d P 2 ( 4> , B ) /dy , 


where <3? and 9 are any two suitably continuous functions 
and P.j ( 4- 9 t : ) and P 2 ( > , B ) are bilinear concomitants. 

By suitably choosing the boundary conditions so that P and 
P 2 vanish at the boundaries, we state the adjoint problem as 

X + A >1*8 = 0 . (3.23) 

Boundary conditions for the adjoint problem turns out to be 
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the same as (3.20) and the adjoint operators are 

X,* = X = dVdy 4 , 

J/U* = 3/2(1 - y 1 2 ) d 2 / d y 2 - 6y d/dy . (3.24) 


If cU and 'd n are eigenfunctions of (3.19) and 

(3.23) respectively'' and the corresponding eigenvalues are 

A and 7\ , it can he shown that 
m n 


( -A 


m 


A ) (/•::„. , 

n J v ’ m 


e,J = o 


■Yhen m / n, it then follows 



U ) = ( do , j P ' 
u n m ' - ti j 


0 


1 

or 6^ dy = 0 , m / n (3.25) 

We can define a new scalar product with respect to the 
operator X » "by the relation 


/n A i 

. ,b J 


= ( :/C ■ho , y ) 


(3.26a) 


where hb and 4? satisfy (3.2 0) and ore' suff iciently 'd ifferen- 
tiahle. Therefore, 


Equation (3.25) then becomes 

, fi | = 0 , m / n 

[_ tn ’ D nJ 

Eigenfunctions of (3.19) and of its adjoint 


(3 .26b) 


(3.27) 

problem (3.23), 
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therefore , form a hiorthogonal set , where the orthogonality 
relation is given hy (3.27). 

3.3.3 EIGENPUNCTIOU EXPAhSIOl 

We will require the expansion of functions satisfying 
suitable conditions into eigenfunctions of the equation 
(3.19). We shall assume that the required eigenfunction 
expansion converges in some sense. 



The perturbation stream function 4^ is small 
sufficiently downstream, so that the nonlinear interaction 
of the modes in (3.3) can be neglected. The solution of the 
line am problem (3.4) can then be written as 


-A (X 


X ) 

° ^J7) * 


(3.30) 


The constant A in the expansion of V. can be found using 
in • 1 c 

the relation (3.29), 
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A 


m 








where 


(3.31) 


. = '"V at x = x 

10 1 o 

3.4 METHOD OF CALCULATION 


(3.32) 


Let in equation (3.2) he written as 

K 

% = X a (x) O m (y) , (3.33) 

m=1 

where the coefficients a^ depend on x, and 4"^ are eigen- 
functions of equation (3.6). In an IT-term calculation 
summation index in (3.33) runs from one to IT . The range of 
summation of indices will not he henceforth mentioned as it 
would he clear from the context. Substituting (3.33) in- 
equation (3-3) followed by multiplication by 8 , and 

integrating, we get 


! 


3, + 


'A 

m 


a 


9 


AT; c 

p" q mpq 


a a„ 

P q 


where 


(3.34) 


Ampq "\i ^impq^Tm 5 (3*35?.) 

^Impq “ ^ "'m * ' r p T'q - ^'p ^'q ) 5 (3.35b) 


r 


2m 


4* , tr 

_ ■ m ’ m j 


(3.35c) 


The differential equation (3.34) governing the coefficients 
of the eigenfunction in the expansion (3.33) is a set of 
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first order ordinary, quasi-line ar , autonomous differential 

f 

equations. 'The first term on the left In equation (3.34) 
represents the linearized part of the inertia term whereas 
the transverse viscous diffusion is represented hy the 
second term. If the nonlinear interaction of the modes 
given by the term on the right side of (3.34) is negligible, 
a linear equation is obtained which has the solution (3.30). 
Equation (3.34) can also be written in the form 


V ( 6 - V C a ) a = - X a 

L.‘ y tnn n mpq a ' p mm 

p - 1 ' 


(3.36) 


The inertia terms have been collected on the left side .in 
the above form. 6^ is the Kronecker delta. 


3.4.1 III IT 111 CONDIT ION 


The flow calculation is now reduced to solving 
equation (3.36) with appropriate initial conditions. 
Initial values of a m can be obtained in a. similar way from 
(3.33). At x = 0 

Mho = I a A 0) X (y) • 

Therefore 


a 


mO 


a m (°) - ( w 10 


. IV 


) / I 


2m 


(3.37) 


If 'm-q prescribed in equation (2.6) is piecewise twice 
differentiable, so is V-jo* 


Equation (3-37) then leads to 
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a mO = 


( ,\F 

( r 0 


) / x 


m 


2m 


3.4.2 LINEAR SUB-RANGE 


(3.38) 


As x goes to infinity, coefficients go to zero. 
However, the computation need not he carried "beyond a certain 
stage when the coefficients are small, in which case the 
solution of the linear problem can be utilised. Then 

a = a (x ) exp j - /• (x - x ) ! , (3.39) 

m m o i_ m o j 

where x is the streamwise coordinate beyond which the non- 
o 

linear term in (3.34) are negligible. 

3.4.3 SYMMETRIC PROBLEM 

When dealing with a symmetric problem, only even eigen- 
functions need to be considered. Index m in (3.33) then 
runs over the set of IT even eigenfunctions. The range of 
integration with respect to y in any integral such as (3.22) 
can therefore be restricted to the upper half of the 
channel, 0 e. 1 

3.4.4 ONE -TERM SO LUT ION 

When only one term is used in the expansion (3.33) the 
resulting differential equation (3.36) can be integrated to 
give 

ln(a.,/a 10 ) - P 111 (a 1 - a 1Q ) = - ?yc , (3.40) 

where a . 0 = a at x = 0. As C„ = 0 for a 1-term 
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asymmetric calculation, we obtain 

- A.x 

a 1 = a 1 o e 1 

3.5 ON EQUATION (3.36) 

Equation (3.36) obtained for the calculation of flow 
development can be written in the form 

C a = - A a , ( 3 . 41 ) 

where a is a vector whose elements are a^ , ag , .... , /I is 
the eigenvalue matrix, and C may be termed as an inertia 
matrix. C and A axe N x K matrices and a is a E vector in 
a calculation based on N terms. 

Equation (3.41) is a set of algebraic equations for 
a , which can be solved using any standard method. A 
typical integration scheme would then be as follows. At a 
given x, the derivative a is obtained from the given a. 
Integration to the next point in x is then carried out using, 
say, the Euler method. 

The derivative a ! can be calculated provided C can 
be inverted, which is possible, only if the determinant of 
C is non-zero. The integration procedure, therefore, 
breaks down for the chosen value of H in the event the 
determinant of C becomes zero at any x. The singularity of 
the resulting equation (3-41)? which arises when the 
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determinant of C is zero, is a general feature of the 
equations of the above tjrpe . Such a difficulty can occur 
when a moment method or an integral relations method is used 
to reduce partial differential equations to ordinary differ- 
ential equations similar: to (3.41). We take up a discussion 
of the phase space to examine further the nature of the 
singularity. 

3.5.1 PHASE SPACE 

An N-term solution can be represented by .an integral 
curve in an N-d imensional phase space, whose coordinates 

are a^ , a 2 , a^. The equation governing such integral 

curves is 

da^/l^ = da 2 /D 2 = = da^/ln , (3.42) 

where D is the determinant of the matrix obtained bv 
m 

replaning the mth column of the matrix C by the vector on 

the right side of equation (3.41). As equation (3.36) is 

autonomous, 3D and D, the determinant ox C, are functions of 
m 

a. . furthermore , they are polynomials of degree IT. 

The origin of the phase space corresponds to the 
fully developed flow and therefore a 1 meaningful' solution- 
should approach this point. The given initial condition 
determines the integral curve in the phase space. The 
evolution of a with respect to x which can be regarded as 

m w 
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a parameter of the integral curve a (x) > can "b e obtained by 
solving one of the equations 

a m = * (3-43) 

The condition D = 0 v;ould normally represent a sur- 
face in the phase space. If an integral curve crosses the 
surface D = 0 (Figure 7), the direction of increasing x 
appears either 'sink' like ur ’source' like and is physi- 
cally unacceptable for describing a real flow. However, an 
integral curve could cross the surface I) = 0 at certain 
singular points of the equation (3.42), if there are any. 

The phase space of symmetric calculation for N = 2 shows 
such a possibility, and is considered here as an illustra- 
t ion. 

3.5.2 PHASE P1AHE 

The axes of the phase plane are a^ and a 2 , and the 
differential equation governing the integral curves is 

da^/da 2 = D.,/D 2 . (3.44) 

Explicit forms of , D 2 > and I) are given in Appendix 2. 
Conic sections that are associated with the quadratic 
expressions, and D 2 , are ellipses and with E is a 
hyperbola. Figure 8 shows the ellipses 1)^=0 and F 2 = 0. 
Evidently, when an integral curve crosses the ellipse 
= 0 or D 2 = 0, except at their intersection, it is 



41 


normal to the or the a 2 axis which is shown Toy vertical 
and horizontal arrows in Figure 8. 

Both D-] = 0 and = 0 pass through the origin 0 and 

are normal to a„ and a 0 axes at 0. 'therefore , 0 is a 
12 

singular point of (3.44) > and it can easily be shown that 0 
is a stable two-tangent node with a^ axis as the tangent. 

Figure 8 shows another point of intersection of 

D = 0 and 3 =0, which is denoted by P. Equation (3.44) 

1 1 

therefore has another singular point P which is shown to be 
a saddle point in Appendix 3. 

The curve D = 0 is also shown in Figure 8. The curve 
is seen to pass through P. This can also be shown from the 
differential equation. If =0 and D 2 = 0 have a point 
of intersection other than 0, at which B ^ 0, then from 

it 

(3.43) a^ and b , are zero. The left side of (3.36) is 
therefore zero, whereas the right side is not. Therefore, 
we conclude that D = 0 at P. 

Given an initial condition, the integral curve in the 
phase plane can be traced. A typical integral curve is 
shown in Figure 8. If the integral curve crosses D = 0, 
except at the second singular point P, the solution is 
physically unacceptable. However at P, as we shall see 
below, it is possible for an integral curve to cross from 
one side of the curve 3 


= 0 to the other. 
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We consider a schematic of a saddle point singularity 
in Figure 9 S where hollow arrows show the direction of x 
increasing if D were positive everywhere. low let us draw 
a curve D = 0 passing through the saddle point, let D y- 0 
above the curve D = 0, and 1 < 0 below it. The direction of 
x increasing is now shown by full arrows. We see that there 
is only one way of crossing from a given side of the curve 
D = 0 to the other that is physically acceptable. Figure 9 
also shows the ’source' or 'sink' like behaviour of the 
integral curves at points other than F of the curve 1=0. 

A point in the phase space which can be joined to the 
origin by an integral curve which is physically acceptable 
will be called accessible. The set of accessible points 
will be called the accessible domain. 

It is therefore seen that all points of the phase plane 
are not accessible. If the given initial condition for 
equation (3.36) is such that the initial point is accessible s 
then no difficulty is encountered in the integration of 
equation (3.36). 


3.5.3 CONCLUDING KEKARKS 

It Is seen in the previous section that some of the 
points of the phase plane are not accessible due to the 
singularity of the differential equation. The occurttrnoe of 
1 = 0 and the presence of the singular points are the 
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features associated with the equation. They can therefore 
occur, in general, with any number of terms in a calcula- 
tion. lor example, in the case of 1-term symmetric 
calculations, the accessible domain is a, 1/C„ . , m - 2.5. 

i U Ml 

3.6 THE INTEGRA! 103 PROCEDURE 

An integration procedure which proceeds in the phase 
space , is adopted with a view to have a simple control over 
the accuracy and the speed. Usual methods of integration 

C 

with x as the independent variable will be applicable when- 
ever the present method is successful. At each step of the 
integration of the differential equation (3-36), the 
derivatives 


s = E , m = 1 , 2 , 
m m 


(3.45 


\ 

/ 


are determined, let n be such that 

|E n j |E m j , m = 1,2, , IT . (3.46) 

Then a is chosen to form the independent variable for the 
n ■ 

next step of integration. Equation (3.36) can then be 
written as 

da ir/ da n = E n/ E n 5 (3.47a) 

dx/da n = 1/E n . (3.47b) 

A simple predictor -corrector method is designed to 
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solve (3.47a). It maintains a uniform level of accuracy 
throughout the calculation. This is effected "by incorpo- 
rating a method of reducing the step length in a » to 
achieve the prescribed level of accuracy when the deriva- 
tives are large, and on the other hand, increase the step 
length to speed up the computation whenever possible . The 
method is self starting. It is described in Appendix 4 
which also includes a brief description of the computer 
program. 


As discussed in the previous section, equation (3.36) 
is singular at P = 0. Therefore, if during the integration, 
D vanishes at any stage , the integration is terminated. 

The above integration procedure is used upto a point 
where the varieties a m are small enough for the solution of 
the linear problem to be applicable. If ja^j -61 f or all m, 
where (a is a prescribed small number, the integration is 
terminated, and equation (3.3S) is used to extend the 
calculation upto downstream infinity. 

3.7 RESULTS OE CALCULATION 

Tho output of the calculation is primarily a (x). 

One can then determine several flow variables. For instance, 
centreline velocity u , pressure gradient dp/dx and wall- 
vorticity co can be determined by 
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u c = 1.5 + £a m ^(0) , (3.48a) 

•dp/dx = -3 + f. a m =i"’ (1) , (3.48 Td) 

and c.- + = ± 3 - >_a m 4 ? m (+ 1) . (3.48c) 

The values of Uj at the upper and the lower walls are denoted 
hy suffix + and - respectively. Stream function is deter- 
mined using equations (3.2) and (3.33). 

Ee attachment point on the wall is determined "by the 
location where wall-vort icity changes sign. In the case 
where reattachment is on the centreline , the change of sign 
of u Q determines its location. The centre of the eddy is 
determined hy finding the location of the extremum of 
stream function within the flow. This operation also yields 
the value of the recirculation. 

In addition, two quantities which are of general 
interest are considered "below. 

3.7.1 DE1/L 10 PMENT LENGTH 

to 

Development length is analog o us ^t he entry length and 
it can he defined as the streamwise distance after which the 
centreline velocity is within a given fraction of the 


ultimate value. 

for example , 




X 99 “ i 

a c ( x q 9 ) ~ 1 -5| 

/I .5 

II 

o 

b 

—A 

(3.49a) 

X 999 ° ! 

U c (x 999 } " 1 * 5 ! 

/1. 5 

= 0.001 . 

(3.49h) 
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That is, x QO is the value of the st re sunwise coordinate at 

J J 

which centreline velocity u differs from its ultimate value 

0 

1 .5 by If-. Similarly x^^ corresponds to a difference of 
O.lf-. 

If x is a streamwise location after which the 
o 

solution (3.39) is applicable, then the development length 
can be found from 

j V a (x ) exp ■- > (x - x ) . !•’ (0) ' = 0.01 x 1 .5 . 

The contribution of the modes higher than the first even 
mode, can usually be neglected, as they decay more rapidly. 
We then have 

Xq 9 = x q - (1/ A g1 )ln (0.01 x 1.5/|a e1 (x o )j) ,(3.50a) 

where the subscript el denotes the first even mode. 

Similar ly 

x qoo = x n - (V ^p^dn (0.001 X 1.5/la ,(x )i) .(3.50b) 

J J J vJ ti I : C I U 

3.7.2 PRESSURE RECOVERY COEFFICIENT 

The non-dimensional pressure gradient is .-3 in 
Foiseuille flow. We define p (x) as 

dp r /dx = dp/dx + 3 . (3.51) 

Ve take p( 0) = p (0) = 0. The pressure recovery coefficient 
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is then defined "by 

= P r (■- ) • (3.52) 

The asymptote to the curve p(x) is therefore given by 

p = -3x + . (3.53) 

If the solution (3.39) is applicable beyond x = x , v 

o ~ r 

is given by 

Vt t 

P x = P r1 + 2.1 a m^ x o ^ A m 5 (3.54) 

vi her a 

F r1 = p r (x o ) = p(x o } + 3x o * (3-55) 



CHAPTER 4 

ENTRY -PLOW HI A C HAMEL s A TEST CASE 


4.1 INTRODUCTION 

We consider the entry flow in a straight channel as 
a; first application of the calculation method developed in 
the previous chapter. 

Although the problem is technically important, the 
attention it has attracted is considerably more than what 
might he expected . Its simple geometry provides a suitable 
testing ground for theories and methods. The present 
application is aimed at establishing the method by comparing 
its results with those obtained earlier by other methods. 

The comparison provides an indication not only of the 
accuracy of the method of integration but also the applica- 
bility of the limit equation and the initial condition 
hypo the sis . 

There are a large number of reported investigations 
on this problem. We, however, consider mainly those which 
are required for the comparison of the results. The first 
anal3 r tical work based on suitable approximations is due to 
Cchlichting (1934* 1968). The flow near the inlet is 
assumed to consist of the wall boundary-layers together with 
an inviscid uniform central region which accelerates to 
satisfy the continuity requirement. A 


series method is 
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used to construct the "bo undary- layer solution in the initial 
portion. The solution for large downstream distances is 
obtained by perturbing the Poiseuille flow, ihe two 
solutions are patched at an intermediate location, this 
solution has been later improved by including additional 
terms in the perturbation series by Collins & Schowalter 
(1962, 1963) and by So id t & Cess (1962). Bodoia & Osterle 
(1961) have obtained the solution of BLE using a finite- 
difference scheme. However, a large truncation error is 
unavoidable , c und further, the scheme is unstable in the 
presence of reversed flow. 

Wang & Longwell (1964) have considered two types of 
entry conditions in their solution of complete USE obtained 
by a f inite-d iference method. The flow is assumed to be 
uniform at the channel entry in the first case and at the 
upstream infinity in the second . The solution for the first 
case is also obtained by C-illis & Brandt H 964 ) and by 
Morihara (1972) . In off-centre maximum of streamwise velo- 
city occurs after a short distance 'downstream of the inlet 
section, and it continues till it is eroded by the thicken- 
ing of the bound ary-layers . This apparently paradoxical 
situation has been clarified in an important paper by 
Van Dyke (1970), who also showed the singular perturbation 
nature of the entry flow problem. The problems considered 
by Schlichting (1934) and by Bodoia & Osterle (1961) are 
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essentially one-term downstream expansions. The first type 
of entry condition appears restricted and introduces frac- 
tional powers of Reynolds number in the expansion. 

4.2 THE INITIAL CONDITION 

The channel consists of two plane parallel plates 
located at y = + 1, x yO (Figure 10). As the flow is 
symmetrical in y, only the upper half of the channel is 
considered . 

The flow calculation as formulated in Chapter 3 is 
reduced to solwing equation (3.36) with an appropriate 
initial condition. At the inlet, flow is taken to be 

u = 1 , 0 kO 7 kl 1 , (4.1) 

= o , y = 1 . 



The number of terms N to be used in a calculation is 
to be chosen in the first place. The ability to predict 
to a good approximation with a few terms is a great advantage. 
We intend to use rather small number of terms in our 
calculation. A 1-term solution is given by (3.40). 
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Two, three and five terms have been used in the calculation 
whose results we discuss in the next section. Calculation 
with more terms is, in principle, straightforward. 


4.3 RESULTS 

The inlet velocity profile is approximated as shown 
in Figure 10. The given uniform velocity profile is approxi- 
mated by functions which vanish at the wall. The region of 
significant error tends to get localised near the wall with 
increasing number of terms. This evidence can be regarded 
as a computational indication of the convergence of the 
approximating functions, although the non-uniformity in con- 
vergence near: the wall is to be expected. A similar situa- 
tion arises when a unit pulse is approximated by trigo- 
nometric functions. To balance the reduced mass flow near 
the wall, velocity elsewhere has to be greater than unity. 
Table 2 gives two indicators of the approximation, namely 

the centreline velocity u at x = 0 and "the associated 

co 

momentum flux M . The momentum flux is related to the inte- 
gral of the square of the error in the streamvise velocity, 

1 

\ (u - I) 2 by = ] u 2 dy - 1 . 

o’ 0 

It is seen from Table 2 that both the quantities approach 
the given value with increasing number of terms . 

I.i . uFUR 

cent .i- uesAgr 

„ M 4 51152. 

\a:. So- 
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IT 

2 

3 

5 

Given 

u 

CO 

1 .206 

1 .122 

1 .090 

1.0 

M 

0 

1 .088 

1 .071 

1 .051 

1 .0 


TABLE 2 RE PEE SENT AT ION OF THE' INITIAL FRO FILE' 


Figure 11 shows the calculated centreline velocity. 

As mentioned shove, u is initially different from unity. 

A convergent trend is clearly seen from the figure. In 
fact , all the three calculations give approximately the same 
u (x), except for a small region near the inlet. 

Pressure grsdient is shown in Figure 12. The three 
calculations are markedly different near the inlet, hut 
slightly away from it the three curves merge into a single 
curve end they approach a common asymptotic value, j dp/dx- 
is expected to he large near the inlet. With increasing 
number of terms, values of j dp/dx \ predicted are indeed 
larger near the inlet, and show a sharper trend. The area 
under the dp/dx curve gives the pressure said is shown in 
Figure 12. 

The results of Morihara (1972) for E = 1000 are also 
shown in Figure 12. The streamwise pressure gradient and 
as well the pressure are functions of both x and y coordi- 
nates in the above calculation. However, they are almost 
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independent of the transverse coordinate for large x. The 
points in Figure 12 are those corresponding to y = 0. The 
present results compare favourably with the pressure gra- 
dient obtained by Morihara. Rapid decay of ; dp/dxj near the 
inlet is approximately represented by the present calcula- 
tion. The pressure distribution of Morihara is in good 
agreement with the 5-term calculation. 

The asymptote to the curve- p(x) is given by equation 
(3.53), where P is the pressure recovery coefficient and 
is the value of (p + 3x) as x-^-oo . represents the 

value of ( p + 3x) at the point upto which the calculation is 
carried (see Equation (3*55) ). F is obtained using 
Equation (3-54). Table 3a contains the values of F ^ and 
P . The difference between the values of any two calcula- 
tions is the area between two respective dp/dx curves. With 
the increasing number of terms, the prediction of dp/dx is 
expected to converge near the inlet with an associated 
reduction of the differences in the numerical values of 
the pressure recovery coefficient. 


E 

"rl 

-P 

r 

2 

0.2187 

0.2437 

w 

0.2645 

0.2913 

5 

0.3245 

0.3505 


TAB IE 3 a 


PRESSURE RECOVERY COEIT ICIEST 
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In Table 3b the pressure recovery coefficient is 
compared with that given by other investigators. Collins & 
Schowalter (1962, 1963) who had earlier obtained a value 

0.338 have later reported 0.37. The calculation of 
Morihara (1972) at E = 1000, gives a value of 0.365. In 
the first case of Wang & Longwell (1964) there is a 


Eeference 

-P 

r 

Schlichting (1934) 

0.301 

Bodoia & Osterle (1961) 

0.338 

Collins &, Schowalter (1962) 

O 

C>1 

VKl 

CO 

" (1963) 

0.37 

Van g & Longwell (1964)" 1 

0.370 2 3 4 


0.352"' 

Morihara (1972) 

0.365^’ 

Pr e se nt 5 -t e r m 

0.351 


TABLE 3b 0 OMULS I SOW OE PEES SUES iiBCOVEEY COEFFICIENTS 


1. Values refer to case II, E=75> and are averages. See P.55 

2. Includes upstream contribution. 

3. Excluding upstream contribution. 

4. For E = 1000. 
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significant variation of pressure across a far -downstream 
section, which is attributed by the authors to the dis- 
continuity at the inlet. Also, an error has been committed 
in the application of the boundary condition an vorticity 
is assumed to be zero at the inlet. In the second case, 
pressure is more or less constant across any section suffi- 
ciently away from the inlet. Three values of pressure 
given at y = 0.1, 0.5 and 0.9 are averaged. Also, these 
values are relative to pressure at x = - CO. ■ Therefore P Vi 
includes a contribution from the upstream of the inlet. 
Thus, the value for excess pressure drop obtained with five 
terms is somewhere in the middle of the range of the values 
given by other workers. We therefore , find that given 
by the present calculation agrees quite well with that 
obtained from calculations based on NSE for large E. 

The wall-vorticity on = - u at the upper wall is 

+ y 

shown in Figure 13. It is also proportional to skin 
friction. Ee suits due to Idorihara for E = 1000 are also 
shown in the figure. The wall-vorticity is infinite at the 
inlet due to the form of the assumed initial velocity 
profile. The wall-vorticity compares well with that 
obtained by Morihara except for a small distance near the 
inlet . 

Velocity profiles at various x based on the 5-term 
calculation are shown in Figure 14. The y show a small 
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waviness near the inlet due to the representation of the 

given velocity profile. The waviness rapidly smoot he ns out 

and the velocity profile appears like a growing wall 

houndary-layer together with a uniform central core, an 

assumed hy Schlichting (1934), (figure 14, see middle row). 

This comparison can he further checked by examining the 

2 

growth of (u - 1) for this intermediate range (,004<x< 
0.02), (Figure 15). The clear linear trend (cf . equation 
(9.69) of Schlichting, 1968) shows that the method brings 
out the qualitative features of the flow. The constant of 
proportionality obtained from the graph is 1 .61 in compari- 
son with the Blasius constant 1.72. The bottom row of 
Figure 14 shows how the velocity prof iles evolve into 
Poise uille flow at distances away from the inlet. 

Development length x ro as defined in section 3.7.1 

yy 

is given in Table 4 along with those obtained by other 
workers. The present values appear to be slightly less 
than those reported by others except that of Schlichting. 
Development length based on the 5-term calculation agrees 
with the value of Morihara (1972) within 1 °jo. 

Figure 16 compares the 5 -term values of u^ with 
those of Bodoia & Osterle (1961), and of Dang and Iongwell 
(1964) for E = 75, and those of Morihara (1972) for E = 1000. 
(Some of the points of Dang £ Longwell shown in Figure 16 
are taken from Figure 3 of V&i Dyke, 1970). The results 
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Reference 

X 99 

Present 2 -term 

0.1642 

3-term 

0.1662 

5-term 

0.1679 

Schlichting (1934) 

0.160 

Hwang & Pan (1964) 

0.1688 

Morihara (1972) 

0.1692 

Bodoia & Osterele (1961) 

0.176 

Gillis & Brandt (1964) 

0.1768 

Roidt & Cess (1962) 

0.1816 

TABLE 4 DEVELOPMENT 

LENGTH 


of Wang and Longwell and of Moxihara agree well with that of 
Bodoia and Osterle near the inlet. However, away from the 
inlet they agree more closely with the present calculation. 
As mentioned earlier, the difference near the inlet could 
he due to the very small number of terms. Possibly 
agreement there could he improved by including more terms. 
Incidentally the limit process used to obtain the governing 
equation corresponds to the downstream limit. Therefore 
the present calculation could be at slight variance, near the 
inlet, with the results obtained for some large value of R. 

A similar agrument also holds for the calculation of 
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Bodoia & Osterle. But however one finds that the agreement 
of their calculation with that based on BSE for large P. is 
better near the inlet than away from it. This difference 
could be due to the accumulation of error as pointed out 
earlier . 

It should be noted that the criterion of determining 
the initial condition leads to a good agreement slightly 
away from the inlet, although the predictions are not as 
accurate near the inlet. That is to say in some integral 
sense the initial condition in the present calculation 
represents the flow well enough for even a few terms to give 
quite accurate predictions slightly away from the inlet. 
Inclusion of more terms may hardly affect the calculation 
away from the inlet (figure 11), however, we expect 
improvement in predictions near the inlet. 

As the basic equations in the present work are the 
same as those of Bodoia & Osterle (1961), comparison with 
their work can be regarded as a check on the accuracy of the 
present method of integration. Also it might be argued that 
a systematic difference in two curves shown in figure 16 
might be due to a. shift in the initial condition. The 
forward-difference formulation of Bodoia t Osterle could 
lead to significant accumulation of error. We therefore 
regard the solution of the complete equations at large 
Reynolds number to be a reliable guide in establishing our 



CHAPTER 5 

APPLICATION TO SEPARATED PLOW PROBLEMS 


5.1 INTRODUCTION 

Having tested the method of calculation for the entry 
flow problem in the previous chapter, we give results in this 
chapter for separated flow problems where many earlier method 
of calculation would not be normally applicable. 

Plow problems in three types of geometrical configu- 
rations are considered (Figure 1). The first type (Figure 
la) is the flow in a channel having symmetric sudden 
expansion. The second type of problem (Figure 1b) is flow 
in an asjTnmetric sudden expansion, which can also be 
regarded as a backward -fac ing step with a top wall, or a 
base with a long splitter plate in a channel. The third 
type of problem (Figure 1c) deals with s base placed 
symmetrically between two parallel flat walls. 

We further consider two types of entry condition in 
each of these problems. A parabolic velocity profile 
corresponds to a fully developed entry condition. A short 
inlet leads to a uniform velocity profile. It Is possible 
to consider intermediate cases having a non-zero boundary- 
layer thickness a little upstream of, say, the base. It 
is also possible to consider other geometrically different 
problems. However, we restrict ourselves to whet seems to 
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"be a representative sample of the above three types of 
problems with one geometrical parameter h and with, two 
simple but extreme entry conditions. 

The number of terms if in (3.33) has to be determined 
for each calculation. We can increase IT till any further 
increase does not affect the results significantly. Alter- 
natively, we can fix N and find the range of the geometri- 
cal parameter in which results are acceptable. 

Consider now the flow in a symmetric sudden expan- 
sion with h = 0.5 when the entry profile is parabolic. ITS 
solution for values of R upto 160 is given by Macagno and 
Hung (1966). Their figure 7 shows that an asymptotic trend 
has been attained and therefore their results can be used 
for comparison with the present calculation. Morihara (1972) 
has also obtained NS solution for this problem for values 
of R upto 20 (see figure 5). The present calculation is 
started with a few terms and later more terms axe included. 
We shall see that five terms are sufficient to indicate 
convergence . Therefore , the number of terms in subsequent 
calculations is restricted to five. 

At this stage it may be appropriate to briefly 
summarise the important steps in the development of the 
present work. A limit equation is obtained for large R for 
a class of problems which may have recirculating regions. 
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A method of flow calculation is then developed which leads 
to a set of or <3 inary differential equations. The upstream 
condition is used directly to obtain the required initial 
condition. The concern of the present chapter is to test 
the consequences of the basic inputs, that is the calcula- 
tion method in its final form when the flows have recircu- 
lating regions. It is applied to a variety of situations 
and comparison with available results is made wherever 
possible. Also some of the results are believed to be new. 

5.2 CHANNEL WITH A SYMMETRIC SUDDEN EXPANSION 

Consider the flow in a channel having a symmetric 
sudden expansion (Figure la). The expansion ratio, which 
is defined as the ratio of the downstream channel width to 
that upstream, is given by 1/h. 

5.2.1 INITIAL CONDITION 


In the case of the parabolic velocity profile 'at the 
entry, the initial flow is given by: 


u q - 3/ (2h) ( 1 - (y/h) 2 ) 


= 0 


, Opyuh , 

» h h Jyj A 1 - 


(5.1) 


The initial condition a_ n is then obtained to be 


m 


, ! . l' 

a _ = -3 > h h - t- 

mO m m 


/ (s„ -n 3 ) 


y=h 


m 


(5.2 



63 


where 

1 

®»- |< e » ^ ‘*hm ■ ( 5 . 3 ) 

0 . 

When the inlet velocity profile is uniform, that is 

u Q = 1/h , 0 < | y l< h , 

= 0 , h <-y 1 . (5.4) 

the initial condition for a is found to he 

m 

a mO = -e" (h) / (G m h) . (5.5) 

We first discuss the case for which NS solution has 
been obtained by Hung and Macagno (1966). Other cases are 
discussed later. 

5.2.2 RESULTS FOE h = 0.5 WITH PARABOLIC ENTRY PROFILE 

The assumed parabolic entry profile given by (5.1), 
is shown in Figure 17 with its approximation using 2, 3 
and 5 terms. Clearly, the given velocity* profile is better 
approximated with increasing number of terms. 3- and 5-term 
representations are closer to each other. 

Streamlines based on the 3-term calculation are 
shown in Figure 18. Also streamlines obtained by Hung and 
Macagno (1366) for R = 46.6 are replotted with the present 
variables in the same figure. It is seen that the stream- 
line pattern obtanned by the present method is in excellent 
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qualitat ive agreement with that obtained by Hung and 
Macagno . The numerical details of the flow are compared in 
Table 5. Values given in the first row are obtained from 
the graphs given by Hung and Macagno. 



Re attachment 

Centre of 

Re c ire ale t ins 


point 

eddy 

mass flow 



X 



r 

0 

: e-1 

Hung and 

Macagno 0.066 

0 . 01 3 

0 . 052 

Pre sent 

3-term 0.064 

0.014 

0 . 045 

TABLE 

5 SYMMETRIC S UDDER EXPANSION , b = 

= 0.5 

Pi* 4 ,, 

JL lx X S 

favourable comparison 

of the preser 

it calcula- 


tion with that of Hung and Macagno (1966) which is based on 
NSE, supports the following views. Equation (2.5) is 
appropriate an the limit equation for large R for the ’far’ 
region, and the streamwise length scale is of order R. The 
method devised to -solve equation (2.5) is applicable even- 
in the presence of recirculating flow. Also the hypothesis 
used to obtain the initial condition to integrate the 
equation appears to give acceptable results. 

The centreline velocity u is shown in figure 19. 

c 

3-term and 5-term calculations give almost the same u 

c 

except very close to the expansion. Although a very small 
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number of terms are used, the results can be said tc be 
indicative of a convergent trend. Velocity profiles at 
various strearawise locations based on the 3-term calcula- 
tion are shown in figure 20. The velocity profile at 
x = 0.06 is approximately that at the re attachment point. 

figure 21 shows the pressure distribution obtained 
using 2, 3 and 5 terms. The point of re attachment is 

marked on each of the curves. At first the pressure incree 
ses followed by a continual decrease as required by the 
Poiseuille flow. It is seen that the maximum pressure 
occurs after the re attachment . The numerical -values of P 
x 99 arL b x x are given later in Table 8a with those for 
other values of h. 

5.2.3 RESULTS POR OTHER VALUES OP h 

numerical calculation with two types of entry 
profiles, the parabolic and the uniform, are carried out 
for various values of h given in Table £. A sero indi- 
cates the occurrence of the singular it].r and the terminatio 
of the calculation (Sections 3.5 and 3.6). h = 1 corres- 
ponds to the entry flow problem. 2-term calculations 
are also carried out for some intermediate values of h 
in addition to those given in Table 6. Table 7 gives 
some typical initial values of a. 




N h = 

2 

3 

5 

2 1 

3 1 

5 1 


0.95 0.9 0.75 0.5 0.25 0.1 

Parabolic entry profile 

1 11 11 1 

1 11 10 0 

1 11 10 0 

Uniform entry profile 

1 11 11 1 

1 1 1 1 1 0 
1 11 10 0 


TABLE 6 CALCULATIONS PEEPOHMEL POP 
SYMMETRIC SULLEN EXPANSION 


1 D ^ 0 

0 : L = 0 for some x (see Section 3.5). 
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a 10 a 20 a 3C a 40 a 50 

Parabolic entry profile 


0.9 

.1322 

-.0160 

.0330 

-.0073 

* 

o 

.Ls. 

0.5 

1 .1668 

-.0200 

-.2516 

.0228 

.0650 

0.1 

3.8088 

1 .0040 

3.1894 

.9347 

2 .3766 



Uniform entry ; 

profile 


0.9 

-.3751 

.0443 

-.0865 

.0176 

.0348 

0.5 

.6744 

-.1758 

.0059 

.1124 

-.1807 

0.1 

3.7262 

.9442 

2.7992 

.7358 

1 .5681 


TABLE 7 INITIAL CONDITION POE 
SYMMETRIC SUDDEN EXPANSION 




68 


For small values of h the initial flow is like a 
strong jet, and therefore the higher modes are expected to 
he significant . On the other hand , for h close to one 
the details of the flow will not he brought out accurately 
unless higher modes are included in the calculation. For 
h = 0.1, the order of. magnitude of (see Table 7 ) 
suggests that more terms should he considered in the expan- 
sion. On the other hand for h = 0.9, the order of magnitude 
of a q suggests that inclusion of more terms may not he 
necessary. However, it turns out that the small 'hack-step' 
requires higher modes to he included for the calculation 
of flow details like recirculation. In the case of h = 0.5 , 
which is in the middle of the range of h, the inclusion 
of the fourth and the fifth modes affect the calculation 
only slightly. 


For h = 0.5 the assumed uniform entry velocity 
profile (5.4) with its approximations is shown in Figure 22. 


The assumed velocity profile has 


disc ont inuit y 


5. 


The approximation is considerably improved when higher modes 
are included, Another example where the assumed velocity- 
profile is discontinuous is the entry flow in a channel 
considered in the previous chapter. However, there the 


d is cont inuit y is 


at the wall, ’whereas in the present case it 


is located within the interval (0,1) 
approx imat ion is improved 'with the i. 
(see also Figure 10). This suggests 


111 loth cases the 
elusion of higher mode 


that discontinuous 


o 
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velocity profiles can "be expanded into eigenfunctions of 
(3-6). figure 22 also shows u obtained with 2, 3 and 5 
terras. The three calculations show a consistent trend. 

figure 23 shows u for three values of h based on 

c 

5-term calculations with both types of entry conditions, 
for the parabolic entry profile the figure shows a monoto- 
nic decrease of u to the Poiseuille flow value. An 
interesting behaviour is observed for h = 0.75 with uniform 
entry profile, first u decreases and then increases to 
attain the Poiseuille flow value. The incoming flow, due 
to its mixing with the ’dead' fluid, first decelerates. 

later as the flow develops, u increases as u (0) u (a.- ) . 

* c c c 

The above qualitative trend is not shown by the 2-terra 
calculation. We therefore see that the inclusing of higher 
modes in the calculation helps in realising the flow picture 
more accurately and thus a. better quantitative prediction 
can be obtained, for h = 0.9 the ’dip’ in u is however 
absent. Possibly the ’dip’ is small and may require more 
terms in the calculation, figure 24 shows the pressure 
distribution based on 5-term calculations for the above 
values of h. While each of the curve for the parabolic 
entry profile show a positive pressure recovery coefficient 
P r , positive as well as negative P occur for uniform entry- 
profile (see Tables 8a, b) . Points of re attachment [ x r have 
been marked in the figure, following x a maximum 
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h 

I'J 

P r 

X 99 

X 

r 

0.95 

2 

.0445 

.0712 

- 


3 

.0521 

.0784 

- 


5 

.0585 

.0706 

- 

0.9 

2 

.1033 

.1157 

- 


3 

.1199 

. 1 1 48 

- 


5 

.1360 

.1140 

- 

0.75 

2 

.3892 

.1876 

.0030 


3 

.4292 

.1862 

.0080 


5 

.4622 

.1860 

O 

CA 

o 

o 

0.5 

2 

1 .320 

.2839 

.0659 


3 

1.234 

.2779 

• 

o 

Cf\ 

v>] 

00 


5 

1.390 

.2775 

.0656 

0.25 

2 

2.022 

.3594 

.1431 

0.1 

2 

1 .838 

.3715 

.1558 

TAB IE 

8a SUMMARY 0E 
PARABOLIC ENTRY 

RESULTS WITH 
PRO EILE 




h 

N 

P 

r 

X 99 

X 

r 

• 1 

2 

-.2437 

.1642 

- 

(Entry 

flow) 

5 

-.2913 

. 1 662 

- 


5 

-.3505 

.1679 

~ 

0.95 

2 

-.2604 

.1673 

- 


5 

-.3147 

.1693 

- 


> 

-.3699 

.1711 

- 

0.9 

2 

-.2700 

.1696 

- 


3 

-.3179 

.1714 

- 


5 

-.3570 

.1726 

- 

0.75 

2 

-.1590 

.1559 

- 


3 

-.1096 

.1571 

- 


5 

-.0955 

.1673 

.0007 

0.5 

2 

.6397 

.1495 

.0219 


3 

.6889 

.1438 

.0234 


5 

.4970 

.1386 

.0220 

0.25 

2 

2.028 

.3432 

.1265 


3 

2 .285 

.3406 

.1253 

0.1 

2 

1 .884 

. 3708 

.1551 

TABLE 

8t> 

SUMMARY OP 

RESULTS 

WITH 


UNIFORM EHTEY PROFILE 
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occurs s as though there is a region of approximately uni- 
form pressure after the re attachment . A similar trend is 
observed in other calculations to be described later. 

Tables 8a and 8b give the va.lu.es of p , x no and x for 
parabolic and uniform entry profiles respectively. 

5.2.4 PEESSUPiE RECOVERY COEPPICIEiiT 

It would be of interest to see how the calculated 
value of P compares with an approximate integral calcula- 
tion. One such approximation is given here. Plow at a 
downstream section is given by Poiseuille flow, and at x = 0 
it is a known flow. The constant pressure gradient of -3 of 
the fully developed flow and any direct viscous effect are 
disregarded. Pressure P acting on the downstream section 
is therefore due to the change in momentum flux. 

We now calculate P with the above approximation. 

r 

Y'hen the entering flow is parabolic (5.1), momentum balance 
between the sections at the expansion and at some distance 
downstream gives 

P r = 6/5 ( 1 / h - 1) . (5.6) 

When the entering flow is uniform (5*3), we get 

P = 1/h - 6/5 . (5.7) 

r 

Approximate calculations have been earlier based on 
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■uniform -velocity profiles which give 

= 1/h - 1 . (5.8) 

It should he noted that P^ in (5.7) could he positive or 
negative depending on the value of h, whereas (5.8) always 
predicts a positive value for P . 

The approximations on which the above formulae are 
based become increasingly inaccurate as h becomes very smal 
or the expansion ratio becomes verj^ large. They can be 
expected to provide an estimate for moderate expansions. 

The calculated values of P are compared with the 

r 

approximate values given by equations (5.6,7) in Pigure 25 
To avoid overlapping some of the points in Pigure 25 (end 
also in Pigure 26) have been staggered. Por a small 
'back-step' i.e. h close to one, with the uniform entry 
profile , P^ is approximately equal to that for the entry 
flow problem (h .= 1). Por the parabolic entry velocity 
profile P goes to zero as h goes to one . It is seen that 
for the expansion ratio upto about two , equations (5.6, 7), 
shown by solid lines, predict P^ rather well. As h goes 
to zero these equations are singular and therefore cannot 
be expected to be valid for large expansion ratios. 2 -term 
calculation predicts an. approximately uniform levc-l of P^ 
for large expansion ratios. P^ is expected to be reduced 
considerably as the losses due to the expansion will be 
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significant for large expansion ratios. It is suggested 
that for the- range of the expansion ratio greater than, say, 
three, P r may he approximately taken to he a constant and 
equal to two- irrespective of the entry condition. 

5.2.5 BE A'i' 'I AC HMEN‘1 AND DEYEbOfMEKl LENGTHS 

The re attachment points x and the development 

r 

length x^ sre shown in figure 26. for h = 0.9 and 0.95 the 
reatta.chment is not obtained with upto five terms. As 
noted earlier, more terms will he needed for smaller ’hack- 
step 1 , to predict accurately the details such as reattach- 
ment. for the ratio of expansion equal to three, Durst 
et . al. (1974) have calculated the flow using USE. Although 
the upstream condition used in their calculation is derived 
from the experiment, it is believed that the entry velocity 
profile is approxima&ly parabolic. I he point of reattach- 
ment as read from their figure 6 is shown. Although R in 
their calculation is only moderately large, x^ obtained by 
them compares well with the present result. Generally x^ 
with the uniform entry profile is smaller compared to that 
with the parabolic entry profile. It is to be expected that 
the f low 7 reattaches earlier in the former case as the 
mixing layer is stronger than in the latter case. x co 
obviously approaches zero with parabolic entry profile as h 
goes to one. However with the uniform entry profile the above 
limit corresponds to the entry flow problem. 
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5.3 CHANNEL WITH AH ASYMMETRIC SUDDEN EXPANSION 

The flow in a one-sided expansion of a channel shown 
in figure 1b is now considered . The primary aim of this 
application is to demonstrate the use of antisymmetric 
e igenfunct ions . 

3.3.1 INITIAL CONDITION 


The downstream channel width is two. The step 

height is denoted "by h. When the entry velocity profile is 

parabolic the initial condition for a becomes 

m 


m 


a m0 


1.5 x (-1) 
(1 - h/2) 3 G 


m 


(1 - h/2) © — 0 j 

m mj y = 1 - h , 


(5.9) 


o £ h £ 1 


where G is given by (5.3). In the case of a uniform entry 


m 


profile, the initial condition for a is obtained to be 


ao) + (-D 


a m0 “ 


m 


’ m 


(1 - h) ) 


(2 - h) G. 


, 0 <£h £ 1 . (5.10) 


m 


Equations (5.9) and (5.10) can be modified when 1 h £ 2 


5.3.2 RESULTS 

Table 9 gives various values of h for which the 
calculations are carried out. With parabolic entry profile, 
the given flow for small h differs only slightly from the 
fully developed flow (see Table 10). The small perturba- 
tion results given in Section 3.3.4 can be used, for h = 0, 
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N 

h = 0.05 

0.1 0.25 0.5 

0.75 

1 .0 



Parabolic 

entry profile 



2 

1 

1 

1 1 

0 

1 

3 

1 

1 

1 1 

1 

1 

5 

1 

1 

1 1 

1 

0 



Uniform 

entry profile 



2 

1 

1 

1 1 

1 

1 

3 

1 

1 

1 1 

1 

0 

5 

1 

1 

1 1 

1 

0 


I 1 ABIE 9 CALCULATIONS PERFORMED 
ASYMMETRIC SUDDEN EXPANSION 

FOR 



1,0: refer 

to Table 

6. 
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h 

a 10 

a 20 

a 30 

P 

O 

k 

i 

a 50 



Parabolic 

entry profile 



1 .0 

-.8595 

-1 .0544 

-.0256 

.0727 

-.0595 

0.5 

-.3877 

.3129 

-.1066 

. 0080 

.0035 

0.1 

-.0108 

.0591 

-.0067 

.0071 

-.0069 



Uniform e 

ntrv profile 



1 .0 

-.6062 

-.3340 

.3203 

. 0401 

-.3017 

0.5 

.1670 

.0021 

-.1056 

.0319 

-.1524 

0.1 

.0659 

-.3535 

.0404 

.0421 

.0405 


TABLE 10 INITIAL CONDITIONS BOR 
ASYMMETRIC SUDDEN EXPANSION 
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the problem is trivial. In the case of uniform entry 
profile, h = 0 corresponds to the entry flow problem. 

Initial conditions obtained for some typical values of h 
are given in Table 10. 

The given parabolic entry profile with its approxi- 
mations is shown in Figure 27 for two values of h. Pith 
increasing number of terms the approximation is seen to 
improve. In the figure instead of showing each of the 
three approximations (2, 3 and 5 terms), only two have beer- 
shown. The omitted, one lies close to the next higher 
approximation. 

When the entry profile is uniform, its approximations 
are shown in Figure 28. Even 5 terms provide a poor appro- 
ximation in this case. In fact for some values of h it 
appears that the total flow is less in the upper half of the 
channel than in the lower half, whereas the opposite is 
true for the given profile. However, a convergence of the 
oscillatory type cannot be ruled out. Here, the numerical 
results for only the parabolic entry profile are discussed. 

The asymmetry of the flow is shown in a simple way 
by considering the skin friction at the two walls. The 
values at the upper and the lower walls .are denoted by 
suffix + and - respectively. The skin friction coefficient 
c f is also related to the wall-vorticity 60 by 
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R = + Uj 


+ 


+ 


\ c f R s.t the two walls are shown in Figure 29 for h=1 . As 
the flow develops the asymmetry of the flow reduces and vani- 
shes ultimately. Compared to the 2— term calculation, the 
3-term calculation shows more asymmetry/ a small recircula- 
ting region is obtained on the lower wall with three terms. 


The point of reattachment is taken to be the point where the 
skinfriction curve meets the x axis going upward (see also 
Figure 30). 


Velocity profiles at various streamwise locations 
based on the 3-term calculation are shown in Figure 30. The 
gradual evolution of the flow into a Poiseuille flow is seen 
from the figure. Pressure distribution is shown in Figure 31. 
Predictions by the 2-term and the 3-term calculations are 
qualitatively similar. At the beginning the pressure drop is 
large followed by a region over which the pressure is appro- 
ximately uniform, and then the pressure decreases further to 
its asymptotic trend. It is interesting to note that the re- 
attachment as predicted by the 3-term calculation marked by x 
in the figure , is followed by a region of approximately uni- 
form pressure. Pressure distribution is shewn in Figure 32 
for some typical values of h. Numerical results are summari- 
sed in Table 11 which gives P,, , x 0 ^ and x £ for various valued 
of h. 

* It should be noted that the eigenvalues are numbered 
in the increasing order, and the first and the third 
e ige nf unc t ions ar e ant 1 s ymme trie. 
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h 

N 


1 

| 

) 

W 

X 

X 

0.05 

2 

.0147 

• 

o 

00 

- 


3 

.0147 

.0318 

- 


5 

.0209 

.0318 

- 

0.1 

2 

.0333 

.0751 

- 


3 

T“ 

C\i 

K\ 

O 

• 

.0732 

- 


5 

.0455 

.0726 

- 

0.25 

2 

.1078 

.1428 

- 


3 

.1017 

.1361 

- 


5 

.1388 

.1319 

- 

0.5 

2 

.3114 

.2101 

- 


3 

.2553 

.1892 

.0094 


r~ 

j 

.3437 

.1883 

.0145 

0.75 

3 

.4104 

.2463 

.0409 


5 

.4252 

.2365 

.0421 

1 .0 

2 

.1617 

.2857 

- 


3 

.1016 

.2799 

.0617 


TABLE 


11 ASYMMETRIC STJEDEN EXPANSION 
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5.4 BASE IN A CHANNEL 


Consider a "base having thickness 2h, placed symmetri- 
cally in a channel as shown in Figure 1c. The flow "behind 
a thick long flat plate with a "bluff end placed in a wind- 
tunnel can "be idealized to lead to the present problem, 
although it is recognised that only a few features of the 
"base flow are modelled in the present problem. 

5.4.I INITIAL CONDITION 


Two types of entry conditions are chosen as before. ' 
For the parabolic velocity profile at the entry, the initial 
condition for a m is obtained to be 


6 


mO 


(1 - h*)G m 


iG m + 2/(1 - h) 0 m 


. (5.11) 


->y 


h 


In the case of the uniform entry profile the initial condi- 
tion for a turns out to be 
m 


1 


a. 


mO 


(i - h)G 
J m 


6 "(h) - e"(i) 

' 'in 


( 5 . 12 ) 


5.4.2 RESULTS 

Table 12 shows the range of the calculations. Table 
13 gives initial conditions for some typical cases. It 
should be noted here that as he- 0 the base shrinks to a 
thin plate. For the parabolic entry profile the problem 



N h = .05 .1 .25 .5 .75 

Parabolic entry rjrofile 
2 1.1111 

3 1110 1 

5 1110 0 


Uniform entry profile 

2 
3 
5 

TABLE 12 CALC IJIAT 1017 S PERFORMED 
FOR THE BASE IN A CHANNEL 


1 1 1 1 1 
1 1111 
10 1 1 0 


1, 0 


; Refer to Table 6. 
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h 

a 10 

a 20 a 30 

a 40 

a 50 



Parabolic entry profiL 



.1 

-1 .4384 

-.0690 -.0769 

.0016 

.0250 

.25 

-1 .7555 

.0738 .2185 

.0296 

-.0377 

.5 

-1 .3457 

.3270 -.5916 

.0205 

.1490 



Uniform entry profile 



.1 

-.7851 

.0604 -.4038 

.0603 

-.2316 

.25 

-1 .3608 

-.0630 -.0646 

.1101 

.1305 

.5 

-1.3424 

.2559 -.1729 

.0737 

.0774 


TABLE 

13 INITIAL CONDITION EOR 
BASE B! A CHANNEL 

THE 
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considered is then like a flow behind s trailing edge of a 
flat plate having finite, equal and opposite vorticity on 
the two sides, When the entering flow is uniform, there 
are two vortex sheers of equal and opposite strength. The 
e igenf unct ions which are used in the approximation have zero 
vorticity at the centreline. Therefore for h = 0, the 
present approximation is non-uniform in a, neighbourhood of 
y = 0. We observe that as h-~* 0, expression (5.12) 
reduces to (4-2) which corresponds to the entry flow problem. 
However for h ^ 0, the approximation is expected to have 
the desired features of the given entry profile such as 
u = 0 at y = 0. 

Figure 33 shows the entry condition as approximated 
by 2 , 3 and 5 terms for h = 0.25 with the parabolic entry 
profile. Table 14 compares the centreline velocity at the 
entry u and the associated momentum flux M given by the 

three approximations with the actual value. Once again we 
see that the approximation is improved with the inclusion ox 


2-term 

3 -term 

5 — ' 1 0 r m 

Actual 

£ 

o 

o 

f 

• 

-A 

00 

.0368 

.0287 

0 

M 1.6139 

1 .5978 

1 .5984 

1 .6 

TABLE 14 REPRESEE 
EOR THE 

TAT ION OF 
BASE IK A 

THE INITIAL 
CHANNEL 

FROFILE 
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higher modes. 

The approx imat ions "bo "uniform entry condition for 
h =0.1 and 0.5 are shown in figure 34. Although any marked 
improvement in the approximation with the inclusion of 
higher modes is not evident for h = 0.5, it is clearly seen 
for h = 0.1. 

'The streamline pattern "based on the 5-term calcu- 
lation is shown in Figure 35 for h = 0.25 with parabolic 
entry profile . A closed "bubble behind the base is obtained . 
Near the base the streamline pattern shows slight waviness 
which is possibly due to the representation of the initial 
condition. This waviness decays rapidly. The development 
of streamwise velocity profiles based on the 5-term calcu- 
lation is shown in Figure 36. The centreline velocity is 
shown in Figure 37. The initially decreasing u c due to 
entrainment is absent in the 2-term calculation, whereas 
3- and 5-term calculations give a better picture of the 
recirculating flow behind the base. Entrainment of the . • . 
'dead* fluid followed by detrainment is seen from the 
figure. Although the size of the recirculating region is 
approximately the same for 3- and 5-term calculations, the 
eddy is stronger and its centre lies closer to the base in 
the 5-term calculation. Table 15 compares the three calcula- 
tions. As the 2-term calculation does not lead to closed 
streamlines, some of the entries in Table 15 are indicated by- . 
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2~term 

3-term 

5-term 

Point of reattachment, x 

’ r 

.0086 

.0086 

.007 

Maximum half -width of the 
recirculating region 

.263 

.205 

.264 

Location of the centre 
of the eddy, 




X 

■ e 

- 

.0042 

.0024 

^e 

- 

.108 

.146 

Total recirculation, - ^ 

-- 

.007 

.021 

Maximum reversed flow 
velocitjr at the centreline 

-.1817 

-.0918 

-.2198 

Development length, Xg 0 

.213 

.206 

.201 

Pressure recovery 

-.492 

-.214 

-.116 


coefficient, p 

’ r 


TABLE 15 BASE IN A CHANNEL (h = 0.25) 
WITH PARABOLIC ENTRY PROEIIE 
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The pressure .distribution is shown in Figure 38. Ths point; 
of reattachment x^ is also marked on the figure . The 
figure suggests that there is a tendency to form a region 
of approximately uniform pressure following the point of 
re attachment . 

Figure 39 and 40 show centreline velocity and pressu 
distributions respectively when the entering flow is uni- 
form. For h = 0.5 and 0.25, u c decreases initially, 
however, no reversed flow is obtained. This evidently is 
due to the error in the approximation of the entry profile. 
It is interesting to note that for h = 0.5 the calculation 
shows a region of 'uniform pressure which in the fewer terms 
calculation is absent. Tables 16a and 16b give P^ , x_g and 
x^ for parabolic and uniform entry profiles respectively. 
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ii 

N 

P 

r 

X 99 

X 

r 

75 

2 

-.4835 

.1846 

- 


3 

-.9255 

.2117 

- 

5 

2 

-.9233 

.2248 

.0216 

25 

2 

-.4921 

.2130 

.0086 


3 

-.2138 

.2057 

.0086 


5 

-.1155 

.2007 

.0070 

1 

2 

-.1621 

.1955 

.0001 



-.2262 

.1947 

.0005 


5 

-.1229 

.1948 

.0006 

05 

2 

-.1596 

.1944 

- 


3 

-.2561 

.1914 

- 


5 

-.1729 

.1916 

- 


TABLE 16a 

BASE IN 

A CHANNEL WITH 



PARABOLIC ENTRY PROEIIE 
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h 

RT 


X 99 

X 

r 

.75 

2 

-.5190 

.1899 

- 


3 

-.9782 

.2157 

- 

.5 

2 

-.8727 

.2253 

- 


3 

-.9321 

.2264 

.0173 


5 

-.3219 

.2244 

- 

.25 

2 

-.1971 

.1970 

- 


3 

-.2360 

.I960 

- 


3 

-.2731 

.1901 

- 

.1 

2 

-.1789 

.1872 

- 


3 

-.4337 

.1822 

- 

.05 

2 

-.2176 

.1795 

- 


3 

-.3854 

.1792 

- 


5 

-.4331 

.1766 

- 


TABLE 16t> BASE IB A CHAMEL WITH 
UNIFORM ENTRY PRO BILE 



EPILOGUE 


EXTERNAL SEPARATEE PLOW 


The internal separated flow provides a suitable 
ground for the testing of limit equations for L, r ^ 1. We 
then find that for '•■■■• 1, two principal limits are possi- 

i 

hie and they correspond to L, r -v' 1 and 1 ^ R. While we are 

"** .a. -K 

familiar with the first limit process which leads to the 
Euler equation, the second limit process which is rather 
new is examined in this thesis. 

The limit equation (2.4) corresponding to the limit 

I rv R is wellknown and is associated with the classical 
x 

hound ar y-layer . We test the applicahilit}’ - of this equation 
to a class of separated flows. Good comparison with solu- 
tions using complete USE supports its validity. 

External separated flows differ in one important 
respect from the internal separated flows. The interaction 
between the outer irrotational flow and the separated 
region is an essential characteristic of the external 
separated flow and poses considerable difficulty in the 
analysis. As it stands, it is difficult to ascertain 
apriori the transverse length scale of the separated wake . 
In fact some drastically different models of the limit flow 
have been proposed. However, If we assume that the trans- 
verse length scale of the separated wake is of order one, 
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some of the conclusions of this work such as the appropri 
ateness of the limit equation (2.4) would he applicable. 

The experimental results given by fcrivos and his 
colleagues (1964* 1965, 1968) for flow past a. circular 
cylinder suggest L 1 and L.'v/ H. The limit equation 

y x 

(2.4) therefore should be applicable. In this limit the 
viscous and the inertial forces are of same order. This 
nature of the flow is different from the other shear 
driven recirculating regions where the limiting flow is 
expected to be governed by the Euler equation. 

The present work may give an insight into one 
element of the external separated flow at large Reynolds 
number whose global picture continues to fascinate many. 
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APPENDIX 1 


SOLUTION .OF THE EIGENVALUE PROBLEM 


1 . The eigenvalue problem obtained in Chapter 3 is 



X ! 3/2(1 - y 2 ) 



( 11 . 1 ) 



0 at y = + 1 



Calculation of odd and even eigenvalues proceed separately. 
The region ox integration is restricted to Os. y 4 1 . The 
boundary condition imposed at y = 0 is therefore as follows 

Odd : X’ = i" ’ =o at y = 0 , (1,1.3) 

!? 

Even ; -p = tp = 0 at y = 0 . (A1 .4) 

For a given A two linearly independent odd solutions 
of (A 1.1 ) can be obtained. One such solution is 



(A1 .5) 


Another solution is obtained on integrating (A1.1) with the 
initial conditions 



A’ l" n' :! r , 
~2 =< h>=T2 =0 at y = 0 



A general solution which is odd can then be written as 
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,(1 ;>\) - 0 . 


(A1 .8) 


Thus the zeros of the function ;X) give? the odd 

eigenvalues. 

Similarly for determining even eigenvalues, we con- 
sider two linearly independent even solutions of ( 1 1 . 1 ) . A 
pair of such solutions is obtained on integrating (A1.1) 
with the initial conditions 

u ? 
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il 
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0 . 
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a non-tr 

ivial linear 
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ci ^ s 
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ie s 
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( : -- 3 
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C.s - -rL - 
' 4 3 

i=,) 

■ 4 
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(A1 . 

11 ) 

1 he 

le 

ft si 

de of (A1 

.11 

) 

is a. function 

of 

X . 

Its 

zero , 


the 

ref 

ore , 

gives the 

even 

eigenvalues. 








The 

values of 

4* Vo 

U i .1 

e 

functions on t 

:ie 

left 

side 

of 


(A1 

.8) 

and 

(A1.11) can 

easily be obtained 

for 

a give 

n 

* 

In 

the 

integration of the 

equation (A1 . 

i) 

the 

fourth 

ord e 

T 


Range -Kutt&-G ills method with s step length 


of 0 . C5 


use a . 


Using the Newt on-Raph son method zeros of the above functions 
are obtained . 


The eigenfunctions are scaled in the 


:c. iicw in 



y = o 

y = 0 
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Odd i = 1 at 

, i 

Even s up = 1 at 


( A 1 .12) 
(A1 .13) 


Solution of the adjoint problem is obtained in a 
similar way. Note that 0 = constant is a solution of the 
adjoint equation (3.23). 'The scaling used is similar to 
(A1.12, 13). 

Odd : 0 = 1 at y = 0 f (A1 .14) 

Even ; tr- = 1 at y = 0 . (A1.1 r 

Eng e invalues are given in Tables A1 end A2 . Each 
time a calculation is to be performed, it is not necessary 
to go through the process of recalculating the eigenvalues- 
They can be given as data in the program and the eigenfunct- 
ions can be generated as solutions of an initial value 
problem. Three of the required four conditions at y = 0 
are given by (A1.3, 12) or (A1.4, 13). The fourth condition 
which can go as data in the program, is also given in 
Tables A1 and A2 . 


m ''hn 

1 U.4335 

2 48.6857 

3 1 04. 572 




( 0 ) 


-5.82126 -8.45962 

5.70202 150.811 

-10.3516 -60.7067 


TABLE A1 


ODD EIGENVALUES 
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n= 

1 

2 

3 

4 

5 

m = 1 

1 

.0003 

.0010 

.0014 

. 0029 

2 

. 0009 

1 

.0058 

.0034 

.0129 

J 

.0009 

.0018 

1 

.0076 

.0086 

4 

.0030 

.0038 

.0193 

1 

.0614 

5 

.0030 

.0036 

. 01 1 6 

.0192 

1 


TABLE A3 

A CHECK 

01 BIORTHOGONALITY 



adjoint pro Diem becomes rather large compared with distance 
over which significant variations occur when higher modes 
are considered. If we decrease the step size, we can expect 
this error to reduce significant ly. However we accept 
these calculated eigenfunctions and efforts to obtain them 
more accurately were not made. Incidentally this exercise 
confirms the development of the biorthogonal set. 


APPENDIX 2 


ON THE PHASE PLANE 


The differential equation (3.36) is written as 


(1 - C 111 a 1 I2 a 2^ a 1 ^ 0 1 2 1 S '1 + C i22 a 2- ): 


, ( A2 • 1 / 


( C 211 a 1 + C 2l2 a 2^ a 1 + ^ ' ~ a 221 a :l C 222 a 2 


)&r>- - \ a 0 . (f 2 .2 } 
a ^2 


tiie above equations can be solved ror and ag ro give 


C 5 ,£i — ZJr^l 9 


vmere D. and 


:e functions of a* end a n £rven "belov? 


(C-i 1 1^221 C 121 C 211' >a 1 


°1 21 C 2 1 1 ^ a 1 + ^° 111 C 222 + C 112 C 221 


°122 C 211 ^ a l a 2 + ^ C 112 C 222 ~ C 1 2 2 °2 1 1 ^ a 2 


l ! i 


222 ^ a 2 + 1 


(A2 .4) 


l u 222 ^2° 121 


' 122^2 


0' 2 ° 1 1 1 A 1 C 212' )a 1 a 2 


(A2.S 


expressions are qudaratic, that is they are of try 
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0 


£| a 1 + /5 2 a 1 a 2 + ^3 a 2 + /Vl + '^5 a 2 + fr 


(.42.7) 


The numerical values of the coefficients, s of B, B^ and 


Dp are given in Table A 4 



jj 

D 1 

jj 0 

r_ 

47 

' 1 

.1796 

-5 .546 

-1 .267 

/■? 

•"2 

-.1793 

74.79 

-22.23 

i-'3 

-3.702 

-611.3 

-203.6 

■<-4 

.6947 

-18.81 

0 

'°5 

<P 'Z 

r\ 

vJ 

-57.55 

'°6 

1 

0 

0 


TABLE A A COEFF IC ISi'II S BI TEE 
SXI3E8SIGF (A2.7) 


7 ^ 


The location of the origin and the 
conic section associated with quadrat 
■ „ and B_ are given in Table A5 . 

I / 


inc 1 inat i on o t 
c expressions 



B 

JJ, 

D,. 

Origin of the 
conic section 

( SU, s &-p) 

-1 .753, 
.3620 

-2 .887, 
-.1766 

2.382 , 
-.2715 

Inclinat ior. in 
degrees 

-1.322 

3-519 

-3.136 

TAB IE A5 J 

DETAILS 01 THE: 

COIIC SSCTI 

:ci S 
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in the transformed coordinates that is the principal 
coordinates of the conic section, D is given as 


D 


0.1817 



3.704 


+ 0.8191 


(A2.S) 


It is seen from (A2.8) that D = 0 represents a hyperbola 
Similarly 1)^ and D 0 are written as 


D 1 = - 3.246 I' 




13.6 9 +27.16 


D 0 = -.6576 t 2 - 204.2 rf + 7.811 


(12.9) 

( 12 . 10 ) 


where 3X0 ^beir respective principal coordinates. 

19 = 0 and D ? = 0 therefore represent ellipses, and are 
shown on Eigure 8. 


2. L short version of the program for the symmetric 
2-term calculation "based on (.42.3) and the coefficients 


.. iven in 2 able JL4> 


be written. 'This is suitable for a 


quick calculation of quantities such as pressure, skin 
friction and centreline velocitv. 



APPENDIX 3 

PEE SECOND SHOuIXii irOH'i'P OP' (3.44-) 

D„ and D„ in equation (3.44) are given by (A2.5, 6) 

i i 

~ue consider the behaviour in the ne igifbo sir ho od of the 
second singular point p. The coordinates of point P are 


■•7309 , a 2p = -.1649 


(A3.1) 


■e cei uif 


>=■ 

-7 


ps 


"IP 


3 .jp So - 


(A3. 2) 


For the analysis in the neighbourhood of P, the 

are considered, be can then 


linear terms in and 


write 


da. 


0 £L 


( jL / ' v * a i I ^ P % 1 * ^ ° ^2 ^ ° a 2 ^ P ^ 2 


s A / o' a> )-> 4- ( o' T) A / j -p 

i ! x i t 4. jr 


(A3. 3) 


yof.pt p, f p 


i 


/ .t T*i / > *5- , /ZhT 

pijpca, y-j ^ t vo- 

i 3 ir > 1 


1 4'S-p ) p g 




(A3. 4a) 


(~ ! Do /raw-,' „ + (E D c /.3 a. ) - f 


2 ; P>2 . 


(A3. 4b) 


Here dash does not denote d if fe rent igt ion with respect to x 
cut with respect to some other parameter.. Inserting the 
numerical values in the above expressions ve get 
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= -24.52 

V5 I 


! 

2 


5 .56 


'>-1 


'1 


- 171.4^ 2 , 

+ 33 .85 ? 2 • 


(13.5 a) 
(13.513) 


The above can also be written as 

.. / 


* 1 ; 

42! 


1-24.52 171 .4 | 

j ; 

I 5 .96 33.891 


! >■— 


'1 , 
2 


(13.6) 


[Therefore trace T and the determinant jQ- cf the matrix on 
the right side of (.43.6) are 


_ Q ^7 

— -/ * J I 9 


-1852 . 


I t -7, 7 % 

V -ti s . * i 


Is .11. <0, P is a saddle point. The eigenvalues ox the 
above matrix equation are found to be 47. S3 and -38.56. 
The principal directions then are, 


'/’I _ 

■11 - 2.56,: 

'?2 = 

0.082 5 1 + 

In t he coord ina t e 

system ( ; 7 

written as 


t 

O = 1 
■'l 

-38.56 i/ . 

/ 1 

_ ! 

7 2 = 

47.93 7 2 


’2 J 
*2 
n ■ 


(13.8a) 

(13.8b) 


(13. 9a) 


(13.9b) 



APIEkDIX 4 

OH 'THE COMPUTER PROGRAM 


1. The predictor-corrector method used to Integrate 

equation (3.47a) is as follows. The predictor used is the 
Euler's method , 



( a m / a n ) H 


(A4.1) 


where b is the predicted value of a. at the next point 

IT* 

of integration, and H is the step length in a . The 
corrector used is the Euler's modified method, 

C a - % = H/2 '(a a /a a ) + (\/yj , (A4.2) 


where c^ is the corrected value. The maximum difference 
"between the predicted end the corrected values is given as 


maxn 

m 



(A4.3) 


In general more than one use of the corrector is 
necessary to get e within a prescribed level of accuracy 
ACT. However the repeated use of the corrector is permitted 
only upto a prescribed number of time MAXOOR within which if 
e does not fell below ACY the step length is halved, E = H/2, 
and the procedure is started again from (i,4.1). However in 


halving the step length if H falls below a prescribed mini- 
mum HlOr , the integration is terminated as the roundoff 


error will be significant . 
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At the end of the integration over the step, equation 
(3.47b) is integrated as 

Ax = H/2 (1/a^ + l/\) , U4.4) 

where Ax is the increment in x. Since the error of the 

3 

corrector (A4.2) is cf order H , the step length is doubled , 
E = 2H, if e 4 ACY/10 = ACI10W , at the end of the integra- 
tion over the step. 

A schematic layout of the above predictor-corrector 
method which incorporates the adjustment of the step length, 
is shown on the next page. At the beginning of the inte- 
gration a..„ and a are known. Equation (3-46) is used to 
determine the independent variable a n • The diagram explains 
the integration over one step. 

Rectangular boxes in the diagram show mathematical 
operations and the logical operations are shown bp diamonds. 
Print statements are deleted to keep the diagram simple. 
Arrows show the flow of operation. The mathematical opera- 
tions are stated within boxes either explicitly or inside 
brackets. Few steps of the diagram are explained. The 

first box contains b (A4.1) which means that b are calcula- 

m m 

ted using equation (A4.1). To obtain the derivative b* , a 

m y 

subroutine Yin IPSE is used, which is shown in another box. 
a ! i 0 i! © iiist i .n ing operations are clear from the diagram. 
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2. The sequence of operations in a calculation is as 

follows; 

Information in the form of whether the calculation is 
symmetric or asymmetric, is required to determine the expan- 
sion set. IT is the number of terms in (3.33). In case of 
symmetric calculation IT is the number of even eigenfunctions 
included in the calculation. 

The required eigenfunctions and its adjoint functions 
are generated. Coefficients and Ig^ given by (3.35) 

are calculated. 

The initial condition for a m given by (3.37) is 
c al cul at e d ne xt , 

The integration of the equation (3.36) forms the 
heart of the program, and is carried out in a separate sub- 
routine 5 PEASEST ’ . This requires a subroutine YFRIME5 to 

, T 

solve the set of equations (3-36) for the derivatives a_. 

The predictor-corrector method used for the integration has 
been described above. In addition following features have 
been incorporated . 

i) If the initial condition is such that B <0, EETURIT 
is executed with a message. 

ii) If at any stage of the computation h changes sign, 
EEfUKfi is executed with 


a message 
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iii) If maxm. !a„ . is less than or equal to a prescribed 
m i m ; 

small number £ , EEfUHE is executed (linear sub- 
range, see Section 3.4.2). 


iv) If the computation proceeds longer than a prescribed 
number of steps, EE1UB1T is executed. 

After the successful execution of PHASE! , the 
desired results are printed out. 




a. SYMMETRIC SUDDEN EXPANSION 



b. ASYMMETRIC SUDDEN EXPANSION 



c. BASE IN A CHANNEL 


FIG. I TYPES OF SEPARATED FLOW 
PROBLEMS CONSIDERED 



FIG. 2 ORDER SPACE REPRESENTATION OF 
LARGE REYNOLDS NUMBER LIMITS OF NSE 

(Log SCALE FOR L x ,L y a Y AXES) 



FIG.3 DOMAIN OF VALIDITY FOR 
LIMIT EQUATIONS WHEN L y ~l 

(Log SCALE FOR L x AND V AXES) 
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FIG. 4 SCHEMATIC OF TWO TYPES OF 
INTERNAL SEPARATED FLOWS 
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FIG. 6 SCHEMATIC OF FORWARD AND 
REVERSED FLOW REGIONS 
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FIG. 8 PHASE PLANE 



















FIG. 15 ( UL C - ! > 2 IN AN INTERMEDIATE 

RANGE OF X 
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FIG-17 REPRESENTATION OF THE ENTRY 
ONDITION FOR SYMMETRIC SUDDEN EXPANSION 


FIG. 18 STREAMLINES FOR SYMMETRIC 
SUDDEN EXPANSION 
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FIG-19 CENTRELINE VELOCITY DISTRIBUTION 


FOR SYMMETRIC SUDDEN EXPANSION 
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FIG. 20 FLOW DEVELOPMENT IN SYMMETRIC 

SUDDEN EXPANSION 





FIG. 23 CENTRELINE VELOCITY DISTRIBUTION 
FOR SYMMETRIC SUDDEN EXPANSION 






FIG. 24 PRESSURE DISTRIBUTION 









FIG. 28 REPRESENTATION OF THE 
ENTRY PROFILE 
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FIG. 29 SKIN FRICTION FOR ASYMMETRIC 

SUDDEN EXPANSION 
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FIG. 32 PRESSURE DISTRIBUTION FOR 
ASYMMETRIC SUDDEN EXPANSION 




FIG. 33 REPRESENTATION OF THE ENTRY 
PROFILE FOR THE BASE IN A CHANNEL 
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FIG. 34 REPRESENTATION OF THE 
ENTRY PROFILE 
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